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INTERACTION OF AN ELECTRON BEAM WITH A FULLY IONIZED PLASMA* 


M. A. Allen and G. S. Kino 
Microwave Laboratory, W. W. Hansen Laboratories of Physics, Stanford University, Stanford, California 
(Received January 30, 1961) 


The amplification of microwave signals on an 
electron beam which passes through a plasma 
has been observed by Boyd et al. and Bogdanov 
et al.* and discussed in many theoretical papers. 
However, in the past, because of the nature of 
the mercury vapor discharge in which these ex- 
periments were carried out, it was difficult to 
obtain a close comparison between theory and 
experiment. Preliminary experiments using an 
almost fully ionized thermal cesium plasma have 
been carried out with encouraging results, and 
are reported here. 

A schematic diagram of the apparatus which 
was used is shown in Fig. 1. An electron beam 
passes through a helix on which a microwave 
signal is induced. The modulated electron beam 
then passes through a 2.2-cm length of cesium 
plasma, which is formed from two hot tungsten 


INPUT 


TRANSDUCER 


spiral filaments immersed in cesium vapor. Ions 
are formed at the filaments by resonance ioniza- 
tion; the electrons are produced by thermionic 
emission. Both the cesium plasma and the elec- 
tron beam are confined radially by means of the 
externally applied magnetic field. The micro- 
wave signal on the electron beam is amplified 
after passing through the plasma and the conse- 
quent enhanced modulation is detected on the 
second helix. Finally, the electron beam is col- 
lected on a collector. A water jacket is provided 
around the plasma region to control the temper- 
ature of the cesium vapor, and hence, the vapor 
pressure and rate of ion emission. Half-wave 
rectified ac was used on the heaters, and all 
measurements of interest were carried out dur- 
ing the off-period of the heating cycle. This 
method of producing a plasma has been shown* 
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to yield a plasma characterized by a very high 
percentage of ionization (greater than 10%), with 
temperatures of the ions and electrons at approx- 
imately that of the heater (2300°K). The dominant 
collision process present is that between elec- 
trons and ions. 

By measuring the gain of the system, from the 
input to the output transducers, with the plasma 
on and the plasma off, it was possible to deter- 
mine the gain due to the plasma alone. The 
parameters of the system which could be varied 
were signal frequency, plasma density, and 
beam current. Assuming that for a given plasma 
density the gain was maximum at the plasma 
frequency, then the maximum gain measured 
with approximately 1 ma in the beam varied 
from 10 db at a plasma frequency of 1000 Mc/sec 
to over 20 db at a plasma frequency of 3000 Mc/ 
sec. At higher plasma frequencies, the gain 
was so high that the tube oscillated at that value 
of beam current. 

A theoretical analysis has been performed in 
order to determine the rf gain for a finite elec- 
tron beam in a finite magnetic field passing 
through a plasma. By making the assumption 
of slow waves, it is found that the fields inside 
and outside the beam have transverse propaga- 
tion constants T, and T,, respectively, where 


1 ~(w */") - lw, » -ku o! 
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and where u, is the electron beam velocity. All 
fields are assumed to vary as exp[i(kz - wt)]; 

wp and Wpp are the plasma frequencies of plasma 
and beam, respectively; and w, is the electron 
cyclotron frequency for the applied magnetic 

field. The assumption was made that the plasma 
frequency of the beam is much less than the cyclo- 
tron frequency, and hence, the third term in the 
denominator of Eq. (1) is negligible. However, 
with the normal values of magnetic field used in 
this type of experiment, such an assumption is 
certainly not justified for the plasma frequency 

of the plasma electrons. Thus, it is felt that the 
theoretical results of Bogdanov et al. which rely 
on this assumption are not applicable to our ex- 
periment. The appropriate matching conditions of 
the fields at the edge of the beam were used to de- 
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termine that for w,.<w<wp», T, is real and T, is 
imaginary. It follows that there should be gain in 
this range of frequencies which becomes infinite 
at the plasma frequency. However, because of 
the nature of the Bessel functions involved, for a 
finite cylindrical beam the gain approaches in- 
finity not as (1 - wp?/w*)-” as in the case of the 
infinite beam, but only as In(1- w»’/w*). Con- 
sequently, any effect in the plasma, such as 
nonuniformity in the density, finite temperature, 
or collision losses, causes the measured gain to 
be finite near w= Wp. Because of the logarithmic 
behavior of the gain parameter, the actual nature 
of the loss mechanism need not be known to great 
quantitative accuracy. 

In the range wpy<w<(w_" +w 2)? the propaga- 
tion constant T, is real, with the result that, if 
the beam is finite and small in diameter com- 
pared with the wavelength, then there are values 
of frequency in this range over which the plasma 
presents an inductive impedance to the beam, 
and gain would be expected. The nature of the 
plasma is not known precisely; thus, it can be 
said only that there is gain within this range, 
and that its magnitude as a function of frequency 
does not vary monotonically with frequency but 
has local maxima and minima. For w>(w¢*+w»')”, 
the predicted gain is zero. 

The theoretical curve for gain in the range 
We<w<Wy for a given wy is plotted with an ex- 
perimental curve in Fig. 2. The effect of con- 
duction currents due to cesium collecting on the 
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FIG. 2. Gain due to beam-plasma interaction plotted 
as a function of frequency. 
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walls of the tube prevented an accurate measure- 
ment of beam current being made. The theoreti- 
cal curves are for an assumed value of 3.5 ma. 

It is noted that gain is obtained roughly from w, 
to (wy +w)™. 

Finally, it is noted that there is an important 
conclusion to be drawn from finite-beam theory, 
which in the range w,<w< (wy" +w,')” gives a 
propagation constant T, which is not only real 
but may be small. This conclusion is, that it is 
possible to obtain gain when the fields radiate 
outwards from the beam, and the wavelength in 
the radial direction is comparable to the wave- 
length in free space. Thus, if the plasma ends 
abruptly, it is possible for such fields to radiate 
directly into free space. This mechanism might 
explain how radiation may be emitted from a 
stellar atmosphere. It may also point to possible 
methods for detecting rf radiation from a plasma 
in the laboratory. 


Further experiments are being carried out 
with a modified cesium plasma tube, and it is 
intended to give a full description of our exper- 
imental and theoretical results in a forthcoming 
paper. The authors are indebted to Sylvania for 
providing parts of a traveling-wave tube and 
to Mr. D. L. Masterson for his part in the con- 
struction of the experimental tube. 





*This research was supported by Air Force Cam- 
bridge Research Laboratories, Air Force Research 
Division, Air Research and Development Command. 
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2E. V. Bogdanov, V. A. Kislov, and Z. S. Chernov, 
Radiotekh. Elektron. 5, 229 (1960) (translated by 
Pergamon Press, New York). 

3M. A. Allen, G. S. Kino, and J. D. Lawson, pre- 
sented at the Conference on Plasma Physics, American 
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AN ERROR IN THE THERMAL CONDUCTIVITY FOR A FULLY IONIZED GAS 


R. W. Vaughan-Williams and F. A. Haas 
Atomic Weapons Research Establishment, United Kingdom Atomic Energy Authority, Aldermaston, England 
(Received January 17, 1961) 


The existence of a discrepancy between the 
thermal conductivities for a fully ionized gas 
perpendicular to a magnetic field as calculated 
by Marshall’ and by Rosenbluth and Kaufmann? 
has been widely realized for some time. Be- 
cause of the current interest in transport phe- 
nomena and the fact that the presence of such a 
discrepancy may lead to doubts concerning the 
validity of the underlying principles of these dif- 
ferent calculations, we feel it important to point 
out the source of this error. 

Rosenbluth and Kaufmann compute transport 
coefficients for large wt, where w and 7~ de- 
note gyration and collision frequencies, respec- 
tively. In particular they find that, assuming 
m;>>m,, the dominant contribution to the thermal 
conductivity K perpendicular to the magnetic field 
comes from ion-ion collisions. This result is 
at variance with the work of Marshall, where in 
the limit of large w7, all types of collisions con- 





1.866w"r? + 0.966 


tribute significantly to the thermal conductivity. 
The other transport coefficients are found to be 
in agreement. 

We used the collision integrals listed on page 
81 of Marshall’ to recompute the moments of 
fy, of Rosenbluth and Kaufmann.’ As expected, 
the first moment leads to Eq. (10) of Rosenbluth 
and Kaufmann but the second moment does not 
give their K. The discrepancy is found to arise 
through an error in the last collision integral of 
Marshall.’ This should read 


[w _2San*(w,”), W _2Sayo'(w2”)] 2 = V2M,9(1 - 6), 


using Marshall’s notation. On inserting this in 
Eq. (5.4) on page 41, we found that a,° and a,’ 
were unchanged, but that 
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and 


wTM, 





yl = 1.25 = 


m, 


corresponding to Eqs. (5.11) and (5.12) on pages 
42 and 43, respectively. 
In the limit of high fields, 


kntkT (‘ 866 + 0.566M, =) 





l=) I= 
all =, =1.25 we 
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This agrees with the result obtained by Rosen- 


wt(w*r? +1.9) 


w"*7"M, +9M, +3.394M," +0.32 wt" +6.282u"7" + 0.933 





bluth and Kaufmann assuming as they do that 
mM; > Me. 





'w. Marshall, Atomic Energy Research Establish- 
ment Report AERE-T/R-2419, 1958 (unpublished). 

2M. N. Rosenbluth and A. N. Kaufmann, Phys. Rev. 
109, 1 (1958). 





NUCLEAR PARAMAGNETIC SUSCEPTIBILITY OF LOW-TEMPERATURE He® 
AND THE FLOW PROPERTIES OF SUPERFLUID SYSTEMS OF FERMIONS* 


Ronald M. Rockmore 
Department of Physics, Brandeis University, Waltham, Massachusetts 
(Received January 16, 1961) 


It has recently been suggested’ that a measure- 
ment of the magnetic susceptibility furnishes a 
convenient experiment for the detection of the 
possible low-temperature phase transition of 
liquid He® to a highly correlated phase. The 
theoretical basis for this suggestion was pro- 
vided by the calculation’ of the temperature be- 
havior of this system in the presence of a uni- 
form magnetic field of arbitrary strength and 
the temperature-dependent magnetic susceptibil- 
ity by the method of Bogolyubov et al.” using the 
model of He® due to Brueckner et al.° 

We wish to point out in this note the simple 
interpretation afforded this calculation by the 
two-fluid model*»® of superfluid many-fermion 
systems and to indicate how measurement of the 
magnetic susceptibility together with the moment 
of inertia of a cylindrical container of the liquid 
determines the effective mass for flow in the 
case of an asymmetric gap.*»* 

The usual restriction to “small fields” (up to 
the order of 10° gauss)* permits us to expand 
Eq. (1) of reference 1 to lowest order in Bx, 
where 8 =1/kT and x is the magnetic field multi- 
plied by y, the appropriate magneton for He’, 
so that one has’ 


x=“*BD {1+ cosh[gs2(k)]} 
k 


: exp[82(k) 
=H'2p 2 ee fees (1) 
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We may also express x simply in terms of the 
statistical factors f(k) as* 


x= 1728) =f (k)[1 -F (k)]. (2) 


Further, since the summand f (k)[1 -f (k)] is 
strongly peaked about the Fermi surface [k=k,], 
one may write® 


x™ (u?/k .*)28 Dek'f (k)[1 -F (&)] (3a) 
~3 u(N am /k,*)=3u7L*[M (0) /k pb (3b) 


where Mf'(0) is a linear combination of the effec- 
tive masses for parallel and perpendicular flow’ 
for the case of the asymmetric gap and is the 
usual effective mass for flow’ in the symmetric 
case. N,, is the total number of normal atoms at 
the temperature T.° Finally, on making use of 
the relation 

X 9 =3H"(Nm /k ”), (4) 
where x, is the “rigid”° or normal value of the 
susceptibility, we have 


(x/Xp) = (N/N )=M 'f (0)/nm, (5) 


where n is the density of atoms. 

For purposes of detailed comparison with the 
work of reference 5, we shall cast the sum 
Lit (k)[1 -f (k)] into its corresponding integral 
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form, namely, 


26 Dos (k)[1 -f (k)] -(L°/E PR \" dxF (x) 
k -1 


in 3 1 
=(L°/E ,\Ky+3K,), (8) 
where 
2 2,-1/2 
Fx)=38(" — aB[B” -10(%)1") [2p PE , PE) 
1C(x)| 
(7) 
and 
1K = ’ ° 
K, +2K, M, (0)/nm; (8) 
as in reference 5, the positive semidefinite 
quantities, K, and K,, are, in units of nm, the 
respective effective masses for parallel and 
perpendicular flow only. Thus, 
(x/Xo) = 3K, + Kz (9) 


Since, as has been noted in reference 5, the cor- 
responding ratio for the moment of inertia, J/J,, 
is equal to K, for the same asymmetric case, 

we have 


K, =2[(x/X9) - U/Ip)] - (10) 


Thus the joint measurement of J/I, and x/x, de- 
termines the effective mass for flow in the asym- 
metric case.® 

It has been noted previously’ that the lack of 
good agreement with experiment for very low 
temperatures indicates that the above result can 
only be considered as valid near the critical 


-2 


, 








temperature. However, within the region of its 
validity, a measurement of the magnetic sus- 
ceptibility should provide a sensitive test of the 
superfluid character of the system since, by 
virtue of Eq. (5), it must be considered one of 
the flow properties of superfluid systems. 





*This work was supported in part by the Office of 
Naval Research. 

'L. H. Nosanow and R. Vasudevan, Phys. Rev. 
Letters 6, 1 (1961). 

2N. N. Bogolyubov, D. N. Zubarev, and Yu. A. 
Tserkovnikov, Doklady Akad. Nauk S.S.S.R. 117, 
788 (1957) (translation: Soviet Phys. —Doklady 2, 535 
(1957)). 

3K. A. Brueckner, T. Soda, P. W. Anderson, and 
P. Morel, Phys. Rev. 118, 1442 (1960). 

‘J, Bardeen, Phys. Rev. Letters 1, 399 (1958); 
R. M. Rockmore, Phys. Rev. 118, 1645 (1960). 

5A. E. Glassgold and A. M. Sessler, University of 
California Radiation Laboratory Report UCRL-9223, 
1960 (unpublished). Our notation follows that of this 
reference closely. 

®This is also the practical range of interest since, as 
has been pointed out in reference 1, fields of the order 
of 10° gauss would be needed to observe the critical 
field even near the critical temperature. 

"K. Yosida, Phys. Rev. 110, 769 (1958). 

*The approximation, Decals dol - -f (kK) = cn (RF) 
xD ¢f (1 -F(&)), where c,,(k) is a polynomial in k, 
is quite usual in calculations in the BCS theory [J. Bar- 
deen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 


108, 1175 (1957)]; see, for example, K. Yosida, Progr. 


Theoret. Phys. (Kyoto) 21, 731 (1959). 
*This definition differs slightly from that given in 

reference 4. 

10R. M. Rockmore, Phys. Rev. 120, 1933 (1960). 
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PHONON- PHOTON DOUBLE-QUANTUM TRANSITIONS AS A DETECTOR 
OF MICROWAVE ULTRASONICS 


N. S. Shiren 
General Electric Research Laboratory, Schenectady, New York 
(Received January 25, 1961) 


We have observed double-quantum transitions, 
in which a microwave photon and an ultrasonic 
phonon, at different frequencies, are simulta- 
neously absorbed. Orton, Auzins, and Wertz’ 
have previously reported similar observations 
in which photons at two different frequencies 
were absorbed in a transition between the S_=-1 
and S, =+1 levels of Ni** in MgO at 77°K. Our 
observations were made at 2°K, on the spin sys- 
tems of Ni** and Fe** in MgO. We are able to 
confirm the double-quantum nature of the narrow 
line observed in the spectrum of the latter,? as 
predicted by Orton et al.’ Further results on 
phonon- photon double-quantum transitions, and 
on the interaction of microwave ultrasonics with 
the two aforementioned ions, will be reported 
more fully. However, it is the purpose of this 
Letter to discuss the use of these transitions in 
the detection of ultrasonic radiation at micro- 
wave frequencies, and to report an experimental 
demonstration of the feasibility of the method. 

Detection of microwave ultrasonics by the con- 
comitant saturation of electron spin resonance 
transitions has been suggested,’ and its feasibility 
has been demonstrated under c.w. conditions.* 
As a detector of microsecond pulses, however, 
the saturation method fails because of the long 
spin-lattice relaxation times which are character- 
istic of paramagnetic impurities in solids. 

Pulse detection of microwave ultrasonics is 
presently effected by the inverse of the piezo- 
electric effect used for ultrasonic generation.°;® 
Because of the small power conversion (~10~‘ at 
X band with a cavity Q =10°) the sensitivity is 
low. The principal disadvantage, however, is 
that the output is not always a measure of acous- 
tic intensity, because it is determined by the 
acoustic wave amplitude, integrated over the 
face of the piezoelectric crystal. The detected 
pulse amplitude is consequently extremely sen- 
sitive to any deviation of this face from being a 
constant-phase surface of the wave. Such effects 
lead to nonexponential decays in pulse echo ex- 
periments, and have been discussed in detail by 
Jacobsen.° 

Referring to Fig. 1, the double-quantum detec- 
tion method consists in the measurement of 
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FIG. 1. Energy level scheme showing a possible 


double-quantum transition for a phonon at frequency 
vy and a photon at vy,. 


microwave power, at frequency v,,, which is 
absorbed in double-quantum transitions, at fre- 
quency v,,, when ultrasonic power is applied at 
frequency Vy =%13- YM: Used as a pulse detec- 
tor, the energy level system is allowed to as- 
sume its steady state under the influence of the 
c.w. microwaves. There is then an increased 
microwave absorption concurrent with the ultra- 
sonic pulse. The method is capable of detecting 
very fast pulses because the rise time for ab- 
sorption is of order T,, rather than the much 
larger time, of order T,, associated with partial 
saturation at low power levels. The decay time is 
also ~T, since, in the absence of double-quantum 
transitions, the return of the 1-3 population dif- 
ference to equilibrium is not observed. Since 
the double-quantum transition probability is pro- 
portional to the product of the two energy densi- 
ties, the microwave power absorbed is a meas- 
ure of the acoustic intensity. By using an en- 
ergy level scheme in which v,,>v,,, and there- 
fore, for a given transition probability, Vy? “y» 
conversion “gain” may be realized, in the sense 
that the rf power absorbed is greater than the 
ultrasonic in the ratio vy/vy." 

We have demonstrated double-quantum pulse 
detection using Fe** in MgO. As a detector of 
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ultrasonic energy the ground-state Zeeman tran- 
sitions of this ion are particularly favorable be- 
cause of the large spin-phonon coupling. Com- 
parative measurements in the same crystal show 
that the spin- phonon interaction is ~ 410° larger 
than that for Mn?*, and ~30 larger than for Ni**. 
The effective spin is? S=1 and, as in Ni**,’ the 

S, =+1 levels are shifted by local departures 
from cubic symmetry, leading to an inhomoge- 
neously broadened line, about 400 gauss wide 
(peak to peak derivative), for the AM =1 transi- 
tions. The 4M =2 transition remains sharp, how- 
ever (6 gauss in our crystals). For a fixed value 
of vy it should therefore be possible to find mag- 
netic field values such that equivalent double- 
quantum transitions will occur over a range 

Vy = Vy * 1000 Mc/sec, and we have observed 
them with Vy = y+ 300 Mc/sec. Used as a de- 
tector, the need for tunable receiving and/or 
generating circuits is therefore eliminated. 

The experimental geometry was similar to that 
previously reported.* A 1.75-cm long x-cut 
quartz crystal transducer was bonded to an MgO 
crystal 1.4 cm in length along a cubic axis. Both 
crystals were 3 mm in diameter. The end faces 
of the quartz and the bonding face of the MgO 
were optically polished. The back end of the 
latter was rough cut. The MgO crystal extended 
across the short dimension of a full-wave rec- 
tangular cavity resonator tuned to 9.2 kMc/sec, 
and the bonded face of the crystal lay just inside 
the cavity. The back end extended outside the 
opposite cavity wall. 0.5-yusec ultrasonic pulses 
were generated at the opposite end of the quartz 
by standard techniques.*® Video detection of the 
microwave absorption was effected by a pair of 
MA408A crystals. 

The observed double-quantum pulses are shown 
on the lower trace of Fig. 2, for a 50-watt ultra- 
sonic drive pulse at 9.45 kMc/sec, and 30-mw 
c.w. microwaves at 9.21 kMc/sec. The upper 
trace is the piezoelectrically detected output at 
the generating end of the quartz. The transmit- 
ter pulse is observed first, followed by echoes 
in the quartz. No echoes from the back end of 
the MgO are seen, due to the rough surface and 
strong absorption by Fe** AM =1 transitions. 

The double-quantum pulses appear halfway 
(measured at the start of the rise) between the 
echoes in the upper trace, since they are ob- 
served at the time of transmission through the 
bond. Their rise and decay times are deter- 
mined by the ultrasonic pulse length, i.e., by the 
time taken for all the energy to be within the 













FIG. 2. Lower trace: double-quantum signal after 
crystal video detection, for ultrasonic frequency 240 
Mc/sec higher than the rf. Upper trace: piezoelectric 
signal after superheterodyne detection. The time be- 
tween successive echoes is 6.2 usec. The initial 
wiggles in the double-quantum output, which appear 
simultaneously with the direct magnetron pulse in the 
upper trace, are due to electrical pickup from the 
pulser. The shape of the piezoelectric echoes is due 
to double pulsing and frequency modulation of the 
magnetron. 


cavity. The total observed double-quantum pulse 
width is thus the sum of this rise time and the 
transmission time through the cavity. The two 
smaller pulses following each large one are con- 
secutive reflections from the back and front 
faces of the MgO. They are observed despite 
the loss of a coherent wave front (by reflection 
at the rough surface) which partially destroys 
these echoes in the piezoelectric detector. Also, 
the nonexponential decay of the echo pattern in 
the upper trace, which is caused by a similar 
effect mentioned earlier, is not observed in the 
double-quantum echo pattern which does decay 
exponentially. The background level is believed 
to be an enhanced rf absorption by those spins 
for which vy =V49 and vy=v9g3, due to partial 
saturation of the 2-3 4M =1 transition by the 
ultrasonics. Its decay time, which should be 
~T,, is shortened by the amplifier response. 

It will be noted that the detection sensitivity is 
lower than for the piezoelectric output with super- 
heterodyne post detection. However, it can be 
improved by using higher rf measurement power, 
better video amplification (noise in the lower 
trace is all amplifier noise), and particularly by 
using higher Fe** concentrations. In our crys- 
tals the concentration is ~0.02%, and the effec- 
tive concentration for the narrow double-quantum 
transition is about 1% of this. The crystal length 
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should be smaller than the one we used in order 
to obtain the best ratio of pulse amplitude to 
total energy absorption. The optimum length is 
one just long enough to completely contain the 
ultrasonic pulse so that the output reaches its 
full amplitude. Results on optimization of the 
detector, together with a detailed analysis, will 
be reported at a later date. 

The advent of an intensity (rather than inte- 
grated amplitude) detector will facilitate experi- 
ments which have not heretofore been possible. 
As an example, it will allow observation of in- 
coherent acoustic energy, important in the study 
of boundary scattering. Absolute ultrasonic en- 
ergy measurements are also possible, since the 
microwave power absorbed may be calibrated if 
the transition probability and numbers of spins 
are known. 

It is a pleasure to acknowledge discussions on 


various aspects of this work with Dr. T. G. 
Castner, Dr. E. H. Jacobsen, Dr. E. B. Tucker, 
and Dr. G. D. Watkins. T. G. Kazyaka aided in 
all the experimental work. 
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4E. H. Jacobsen, N. S. Shiren, and E. B. Tucker, 
Phys. Rev. Letters 3, 81 (1959). 

5E. H. Jacobsen, Quantum Electronics (Columbia 
University Press, New York, 1960), pp. 468-482. 

6H. E. Bommel and W. Dransfeld, Phys. Rev. 117, 
1245 (1960). 

'E. H. Jacobsen has pointed out that the inverse, in 
a properly chosen energy level scheme such that 
Vy>> Yy,, allows detection of very high frequency 
phonons by low-frequency rf. 








SURFACE STATES ON CLEAVED SILICON* 


D. R. Palmer, S. Roy Morrison, and C. E. Dauenbaugh 
Honeywell Research Center, Hopkins, Minnesota 
(Received November 11, 1960) 5 


Measurements of surface band structure on 
silicon (111) surfaces cleaved in vacuum (10° to 
10-*° mm) have been made using both the field- 
effect approach’»* and the channel method of Statz 
and his co-workers.* 

In the field-effect measurements the portion 
cleaved off the sample was used as the field 
plate, and measurements were taken with fields 
up to 10* volts/cm. The field-effect measure- 
ments show a slightly n-type surface on both p- 
and n-type bulk material, of field-effect mobility 
of about 50 cm?/v sec. The surface is extremely 
insensitive to oxygen admission, with up to at- 
mospheric oxygen pressure producing negligible 
change in field effect, and no sudden changes of 
conductance with increasing pressure. The lat- 
ter measurement is unreliable over long periods 
of time due to temperature drifts. Because of 
this insensitivity to ambient, the measurement 
cannot distinguish between a strongly n-type 
surface with a high density of states or a surface 
which is almost intrinsic. 

To discriminate between the above techniques, 
the channel method of Statz was employed. Sam- 
ples were made by diffusing phosphorus into 
80-ohm-cm p-type silicon to form npn units and 
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by diffusing boron into 20-ohm-cm n-type silicon 
to form pnp units. The “base” width was 0.04 
inch. By cleaving so as to leave the junctions 
intersecting a cleaved surface, we were able to 
make measurements of the channel conductance, 
if any, caused by the cleaved surface. In neither 
the npn nor the pup case was the conductivity of 
the cleaved surface greater than 10°? .mho/ 
square at 1.5 volts bias, indicating there was no 
channel (inversion layer) present. This can be 
compared, for example, with a conductance of 
about 90 ,.mho/square as measured by the Statz 
technique on a cleaved surface of 5-ohm-cm n- 
type germanium (pnp configuration). 

Thus we conclude that the dominant surface 
states are near the center of the gap with a 
freshly cleaved silicon surface, if we assume 
the surface state structure independent of the 
bulk type. 

This result is in disagreement with work on 
ion-bombarded surfaces, where Law‘ has found 
a p-type accumulation layer on p-type silicon 
after bombardment and annealing and has found 
that the surface becomes n-type after bombard- 
ment with 1000-volt argon ions. He has also 
stated without reference that earlier work showed 
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a p-type surface on 7-type silicon. Dillon and 
Farnsworth® have concluded via work function 
and photoelectric threshold measurements that 
a degenerate p-type layer is found on a silicon 
surface cleaned by ion bombardment and anneal- 


ing. 
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STORED ENERGY RELEASE BELOW 80°K IN DEUTERON-IRRADIATED COPPER* 


A. V. GranatoT and T. G. Nilant 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received December 22, 1960) 


Experiments by Blewitt, Coltman, and Klabunde’* 
on neutron-irradiated copper have shown that the 
energy released during annealing in the vicinity 
of 40°K, designated Stage I, is 1.6 calories/g 
per pohm cm of resistivity which anneals. This 
value is anomalously low and leads to doubts as 
to whether or not interstitial-vacancy mutual 
annihilation can be the recovery mechanism. 
Further, since the annealing observed in deu- 
teron-irradiated copper is in some ways similar 
to that found in neutron damage, the doubt at- 
taches also to the deuteron results. An experi- 
ment was therefore performed to measure the 
energy released after 11-Mev deuteron irradia- 
tion. 

For the measurement, a specially designed 
differential calorimeter was constructed in which 
it was possible to measure all heat leaks as a 
function of temperature. This was required 
since the total anticipated energy release is 
small, and leaks are a significant part of the 
total effect. The magnitude of the energy re- 
leased was limited on the one hand by the small- 
ness of the specimen and on the other hand by 
the necessity to keep the total stored energy re- 
leased per gram per temperature interval con- 
siderably less than the specific heat to avoid 
making the specimen thermally unstable.” The 
size of the specimen was limited by the deuteron 
beam size and range. The above requirements 
were met by using a }-g specimen,® 5x10~ in. 
thick, in deuteron fluxes several times smaller 
than those which have been used in previous 





radiation-damage measurements and a relatively 
fast annealing rate of about 2°K/min. The heat 
leaks were minimized to an extent consistent 
with the requirement of reproducibility by using 
high thermal resistivity Advance-Chromel P 
thermocouple wires and pressures less than 107° 
mm Hg in the experimental volume. 

Departures of the differential system from 
perfect symmetry were taken into account by 
use of a separate background annealing run, in 
which the specimen in an unirradiated state and 
its symmetric dummy were heated to 80°K, and 
subsequently cooled. This run was also used to 
measure the heat leaks and the specific heat of 
the specimen as a function of temperature. The 
important influence of the rapidly varying specif- 
ic heat on the results as well as other aspects 
of the measurement will be discussed in detail 
in a later publication. A continuous recording 
of the output of a differential thermocouple be- 
tween the copper foils was made during the anneal. 
From this recording, calculations could be made 
at temperatures as closely spaced as desired. 

The results are shown in Fig. 1, where the 
energy released per gram per unit temperature 
is plotted in millicalories per gram-°K as a 
function of the annealing absolute temperature. 
Two irradiations and subsequent anneals were 
made: Run I to an integrated flux of 8.25 x10"* 
deuterons/cm? and Run II to 2.89x10** deuterons/ 
cm’. The part of the curves below 27°K is blurred 
by a complication due to experimental difficulties 
in fixing the starting point of the annealing, and 
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FIG. 1. Stored energy 
release in deuteron-irra- 
diated copper. e-—Run I; 
@=8.25 x 10° deuterons/ 
cm’; (E/m)25° - 55°K = 92.8 
millicalories/gram. a- 
Run II; =2.89x10" deu- 
terons/cm’; (E/m)o5 _ 55° 
= 35.3 millicalories/gram. 








should be ignored at this time.* Between 27 and 
55°K there are three peaks. In Run I, the large 
peak has a pronounced high-temperature shoulder 
at about 50°K. Between 55 and 80°K there is a 
more or less structureless additional annealing. 
When the two runs are normalized to the same 
total flux, the three main peaks are in close 
agreement as to size and shape, except for the 
high-temperature side of the large peak. 

In Fig. 1, one finds 92.8 and 35.3 millicalories/ 
g as the total energy released between 25 and 
55°K, for Run I and Run Il, respectively. To 
facilitate computations and comparison with the 
earlier deuteron experiments, these results 
were normalized to a flux of 107 deuterons/cm? 
and a deuteron energy in the foil of 10 Mev. An 
11-Mev deuteron has an average energy in 5 
mils of copper of 7.04 Mev. The average of 
the normalized energy release is then 0.83 cal/g. 
The claimed accuracy has not yet been fixed, 
although it is expected to be of the order of 10% 
or less. To relate this energy release to pre- 
viously determined experimental observations 
of resistivity and volume changes, we use the 
resistance increment given by Cooper, Koehler, 
and Marx,® the fractional recovery found in 
Stage I by Simmons and Balluffi® (less the low- 
est temperature peak) of 54%, and the volume 
change found by Simmons and Balluffi and Vook 
and Wert.’ The resulting experimental ratios 
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normalized to the same deuteron flux and energy 
are [A(E/m)]/Ap=7.1 cal/g wohm cm and (AE/AV) 
= 4.2 ev per atomic volume change. If the energy 


per Frenkel pair (F.P.) is assumed to be 5.0 ev, 
these figures would lead to a resistivity of 2.6 
ohm cm per 1% F.P. and 1.2 atomic volumes 
per F.P. Meechan and Sosin® recently found a 
value of 5.4 for [A(E/m)]/Ap in a measurement 
of the stored energy released at low temperatures 
after electron irradiation. The sensitivity avail- 
able to them was less than that used here, and 
they indicated that their results could well be 
higher. 

When the curves of Fig. 1 are compared with 
the recent isothermal resistivity measurements 
at higher doses of Magnuson, Palmer, and 
Koehler,® one finds a close agreement for Run I 
for the peak structure in Stage I. The peak 
positions in the present work are shifted some- 
what to higher temperatures because of the high- 
er heating rate. The peak shapes are better de- 
fined in the energy measurements because the 
data are essentially continuous. This better 
definition shows structure in the large recovery 
peak at 42°K, present in both runs. Although the 
additional continuous annealing above Stage I 
(above 55°K) is not observed in the resistivity 
measurements, an additional annealing in this 
range is found in the lattice-parameter meas- 
urements of Simmons and Balluffi. 
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When the results are compared with those of 
Corbett, Smith, and Walker’® for electron irra- 
diation, one sees again agreement as to number 
and location of peaks in Stage I, although little 
or no Stage II is observed in electron irradiation. 
Corbett, Smith, and Walker and Koehler, Seitz, 
and Simmons” have interpreted peaks A (below 
20°K), B (~29°K), and C (~34°K) as close inter- 
stitial-vacancy pairs, peak D (~43°K) as cor- 
related recombination of an interstitial with its 
own vacancy, and the concentration-dependent 
peak E (~50°K in Run IJ) as uncorrelated long- 
range interstitial recombination. Our results 
support this interpretation of peaks B, C, and 
E, and, in addition, permit further conclusions 
to be drawn with regard to them and to peak D. 
The close-pair peaks of Corbett, Smith, and 
Walker have half-widths very close (~2°K) to 
the natural half-widths which are to be expected 
for a single discrete activation energy and low 
doses, according to calculations we have made 
which will be presented later. This analysis, 
which takes into account the thermal instability 
mentioned earlier,” shows that the half-width of 
the energy release curve is increased when the 
energy release per degree is a sizable fraction 
of the specific heat. This condition is readily 
obtained at low temperatures and typical deuteron 
damage levels. The half-widths we find for 
deuteron irradiation are somewhat wider than 
those which should occur for a unique activation 
energy (see, for example, peak B). This can 
also be taken as evidence for clusters in deuter- 
on damage since the resulting strain fields ina 
cluster could broaden out a unique activation 
energy somewhat. In addition, one finds that 
peak D is far too wide for a unique activation 
energy but can be interpreted as a superposition 
of many close pairs with activation energies 
closely spaced and approaching a limit which is 
equal to that associated with the freely-migrating 
interstitial of peak £. This conclusion is sup- 
ported by the evidence of structure in peak D. 
Thus the large width of peak D, plus the fact 
that peak D has structure, leads us automatical- 
ly to the close-pair model of radiation damage 


at low temperatures. 
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RESTORATION OF STABILITY IN FERROMAGNETIC RESONANCE 


H. Suhl 
Physics Department, University of California, La Jolla, California 
(Received January 30, 1961) 


As has long been recognized,'~* ferromagnetic 
resonance is subject to instability when the 
radio-frequency signal exceeds a certain thresh- 
old value. In the more common type of instabil- 
ity, pairs of spin waves with equal and opposite 
momenta grow far beyond their normal, thermal 
equilibrium amplitudes at the expense of the uni- 
form precessional motion induced by the signal 
field. In this paper we show how, by rather sim- 
ple means, this threshold may be substantially 
raised. The method is reminiscent (in the time 
domain) of the alternating gradient stabilization 
(in the space domain) of charged beams. Studies 
of large motions of the magnetization, which 
have proved very fruitful in nuclear resonance, 
may now become possible. 

The method we propose is based on the follow- 
ing observations. The natural frequency of the 
spin-wave pair prone to unstable growth stands 
in a very sensitive relation to the signal fre- 
quency, as is obvious from energy conservation 
requirements. If now the natural frequency of 
these spin waves is varied in time, no pair of 
waves will satisfy the necessary frequency rela- 
tion for longer than an instant, and this must in- 
evitably cause an increase in the threshold sig- 
nal. The extent to which this principle is borne 
out in practice is described in another Letter in 
this issue.® 

In standard form, the equations of motion for 
a potentially unstable spin-wave pair of wave 
numbers (k, -k) may be written 


Here wy is the natural frequency of the spin 
wave; af anda@p, respectively, are the ampli- 
tudes of the spin wave and of the uniform preces- 
sion. Py is a coefficient measuring the strength 
of the nonlinear coupling. Ui is a phenomenolog- 
ical loss constant for the spin wave. n is a posi- 
tive integer measuring the order of the nonlinear 
coupling. Only the cases n=1, 2 have been ob- 
served so far. Equations (1) predict exponential 
growth of the spin-wave pair when 
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n 
la o! > ng/ Pe, (2) 


provided the frequency w which drives a, is such 
that 


w= 2we/N. (3) 


Equation (3) determines a surface in k space, 
upon which the right-hand side of Eq. (2) is mini- 
mized, to find which spin wave consistent with 
(3) goes unstable first, and at what minimum 
power. 

We now impose a slow sinusoidal modulation 
of frequency v and amplitude / on the applied 
steady field. v will ordinarily be much less than 
w, OF wy. So long as this is so, the equations of 
motion remain very nearly the same as (1), ex- 
cept that the constants therein become sinusoidal 
functions of the time. For reasons which will 
become clear, the time variation of P; may be 
neglected. The time variation of a, will be im- 
portant only in the case n=2, when the modula- 
tion will tend to sweep the field clear through 
the resonance line, and then only if the modu- 
lating field exceeds the linewidth. In this paper 
we neglect the time variation of a, which is 
easily taken into account, if desired. 

To first order in the modulating fieldh, we 
have, for the time-dependent spin-wave fre- 
quency, 


we’ = Wet (Bw -/ aH )h cosvt. (4) 


A change of variables, 


a= =Ay expti[wt+ (h Owe/ veH)sinvt| - nt} ; 


Ay =, exp(iwt), 
renders Eq. (1) of the form 
a. n ies ees _* 
Ag= PA, exp|i(nw 2we)t 2i(f/ v)sinvt]A » (5) 


together with the “adjoint” equation. Here we 

have defined an additional frequency, f =hdwy/ 8H. 
Let us now confine further discussion to the case 
n=1, for definiteness, and also because the neg- 
lect of the time variation of A,= |a,| is then most 
certainly justified. The condition most favorable 
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to instability is then, provisionally at any rate, 
we=w/2. For spin waves satisfying this condi- 
tion, we have 


Ag=A PeA ¢* exp(-2if/v) sinut 
+00 


=A o£ u Jn (2f/v) exp(-imvt). (6) 
The right-hand side of this equation consists of 
one secular term and a series of nonsecular 
terms. In the crudest approximation only the 
secular term contributes to the instability phe- 
nomenon. In this approximation, then, we retain 
only the term with the secular coefficient, PEdp. 
Instead of the instability criterion (2), we then 
have 


Hence if the modulating field amplitude and the 
modulating frequency are chosen so that the ar- 
gument of the zero-order Bessel function in this 
expression is equal to one of its zeroes, the 
threshold becomes theoretically infinite. Ina 
more refined approximation the threshold is, 
generally, finite. The solution of (5) and its ad- 
joint will have the form 


at. . iret 
-= B 
Ay ° ye vr . 


where the B’s in the summation are constants, 
satisfying the equation 
+00 
~ 2 ° 
B =A) IPE melee Ms Mr)B,» (7) 


where 


+00 


Kiu,m;rx)= DJ. J 


Atirv)(At+inv 
enor pm! Driv rine), 


the argument of the two Bessel functions in this 
expression being the same as above. Since 

Eq. (7) leads to a very complicated infinite de- 
terminant relation for :, we instead discuss a 
completely soluble case: that of square-wave 
modulation. For the first half-period, 7/2, we 
assume that the frequency condition w-=w/2 is 
satisfied exactly; for the second half- period we 
assume a deviation f. Then the equation for A 
[= ag exp(iw7t) | becomes 


(0<t<$T) 
($T7<t<T) 


A=A,PA*, 
A =ifA+A,PA*, 


where the k subscripts have been omitted. 


Equations (8) are solved in the two half-per- 
iods and fitted together at $7 by the requirement 
that A,A* are continuous there. Next, a linear 
combination of A(T) and A*(T) is formed in such 
a way that it is proportional, with a factor of 
proportionality exprT, to the same linear com- 
bination of A(0) andA*(0). Stability is achieved 
if r< n- Some algebra shows that ) is given by 


AT = cosh" ‘[cosh}P,T cosh}2T 
+(P,/Q) sinh}P,T sinhiQ7], 


where Q?=P,?- f? and where P,=A,P. A partic- 
ularly interesting condition arises when the ar- 
gument of the inverse hyperbolic cosine is less 
than unity. Then \ is imaginary and the system 
is stable, no matter how small nj. It is easily 
verified that this condition of absolute stability 

is satisfied only when 2 is imaginary, i.e., when 
the depth of modulation f exceeds A,P. This is, 
of course, to be expected from the general thresh- 
old conditions (see references 1-5) for spin waves 
failing to satisfy the frequency condition. For 
modulation frequencies much smaller than A,P 
the absolute stability criterion is particularly 
easily interpreted. First, for general T, abso- 
lute stability is achieved when 


|cos($2’T - ¢)|<R™}, 
where 
R =[cosh?(3P,T) + (P,?/2”) sinh?($P,T)], 


N27 =f?-P,?, (>0) 
and 
tan $ = (P,/2’) tanh($P,7). 
For large P,T we thus require 
£2'T > p= (n+4)z. 


The interpretation is that for large P,T only the 
growing solution of (8) is of importance at the 
end of the first half-cycle. For f>P,, the solu- 
tion in the second half-cycle is purely oscillatory 
and, by judicious choice of 7, etc., can always be 
made to vanish at T, no matter how large its 
amplitude was at 37. 

All the results in this note apply to instabilities 
due to parallel pumping as well. It is necessary 
only to redefine A, as the pumping field ampli- 
tude itself. 

Stabilization can also be achieved by modu- 
lating the signal frequency, but this will gener- 
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ally be less convenient. 
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SUPPRESSION OF SUBSIDIARY ABSORPTION IN FERRITES BY MODULATION TECHNIQUES 


T. S. Hartwick, E. R. Peressini, and M. T. Weiss 
Components Division, Hughes Aircraft Company, Culver City, California 
(Received January 30, 1961) 


In the preceding paper,’ Suhl has shown theo- 
retically that the threshold power for subsidiary 
absorption in ferrites can be increased by super- 
imposing on the dc magnetic field, Hg., a small 
ac modulation, H,,. Experiments have been 
performed which confirm the above prediction 
and show that subsidiary absorption can indeed 
be suppressed at power levels as high as 10 db 
above the threshold power. 





H,, ARBITRARY UNITS 








i = aes | 


10 
fms Me 





FIG. 1. The minimum modulating field, H,,, re- 
quired to suppress perpendicular pumped subsidiary 
resonance as a function of the modulation frequency 
for various microwave fields. 
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The experiments were performed on an yttrium 
iron garnet (YIG) sphere 0.059 inch in diameter 
with a uniform precession linewidth of 1.2 oe. 
The sphere was placed near the wall of a cavity 
resonant at 9250 Mc/sec. Eddy current losses 
were minimized by using a cavity wall 0.002 mm 
thick and thus an external modulating coil was 
suitable. In Fig. 1 a plot is shown of the relative 
amplitude of H,, required to suppress subsidiary 
resonance for various values of microwave field, 
h. For these data h was oriented perpendicular 


(a) Hy, = 0.00 oe 


(b) Hy, =0.18 0€ 


(c) Hy, =0.30 oe 


(d) H,, =0.36 oe 


FIG. 2. Absorption is recorded on the vertical axis 
as Ha, is traversed through parallel pumped subsidiary 
resonance. Horizontal scale is 20 0oe/cm. The modu- 
lating field is increased from trace (a) through trace 
(d). 





VoLUME 6, NUMBER 4 


PHYSICAL REVIEW LETTERS 


FesBruary 15, 1961 





to Hgc which was set to 2054 oe, the subsidiary 
resonance peak. The results show that the opti- 
mum modulation frequency is about 0.8 Mc/sec 
for a microwave field near the threshold field, 
he. Furthermore, it is evident that H,, is quite 
sensitive to changes in either modulation fre- 
quency or microwave field. These results are 
in general agreement with the Suhl theory al- 
though a detailed comparison cannot as yet be 
made. 

The above Suhl theory can, of course, be ex- 
tended to the parallel pumped case where the 
microwave magnetic field is parallel to the dc 
field. Figure 2 shows the results obtained for 
this parallel pumped condition. In Fig. 2(a) an 
oscilloscope trace shows the subsidiary absorp- 
tion as a function of Hg. with no modulation ap- 
plied. The noisy absorption is due to the relaxa- 
tion oscillations which usually accompany sub- 
sidiary resonance.” As the modulation field 
amplitude is slowly increased, the subsidiary 


absorption gradually becomes suppressed as 
shown in the succeeding traces 2(b), 2(c), and 
2(d). For this microwave field (h/h,=1.26) a 
modulating field of 0.36 oe is required for com- 
plete suppression. The optimum modulating 
frequency is 0.25 Mc/sec as compared to the 
0.8-Mc/sec frequency for the perpendicular 
pumped case of Fig. 1. At power levels of 5 db 
or more above threshold the disappearance of 
the absorption with increasing modulation ampli- 
tude takes place in a discontinuous manner and 
will be described in detail in a forthcoming paper. 
The application of these results as a tool to ex- 
plore the spin-wave spectrum and to various 
microwave devices will also be discussed. 





'H. Suhl, preceding Letter [Phys. Rev. Letters 6, 
174 (1961)]. 
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OSCILLATORY MAGNETORESISTANCE IN n-TYPE PbTe 
C. D. Kuglin, M. R. Ellett, and K. F. Cuff 


Lockheed Research Laboratories, Palo Alto, California 
(Received January 6, 1961) 


Distinct oscillations of the transverse magneto- 
resistance of three oriented single crystals of 
n-type PbTe in dc magnetic fields have been ob- 
served. These results confirm the existence of 
a multivalley band structure having prolate en- 
ergy ellipsoids oriented along the (111) crystal- 
line axes.’»? It is further established that there 
are four ellipsoids, thus locating the energy 
minima at the zone edge. The temperature de- 
pendence of the oscillatory amplitudes yields a 
value for the cyclotron mass m,* and, in addi- 
tion, the transverse and longitudinal effective 
masses are calculated. A recent note,® describ- 
ing the behavior of the transverse magnetoresist- 
ance in pulsed fields in an n-type PbTe crystal, 
has given indications of oscillatory effects. How- 
ever, the measurements, showing a single mini- 
mum, were not interpreted in terms of a possi- 
ble model for the band structure. 

The bulk of our measurements were made on a 
[110|-oriented sample having a carrier concen- 
tration m of 3.45x10'"/cm® and a Hall mobility of 
1.2 x 10° cm?/volt-sec at 1.7°K. With the current 
along the [110] direction, measurements as 


shown in Fig. 1 were obtained with magnetic 
fields up to 26.5 kgauss in the [001] and [110] 
directions at 1.7 and 4.2°K. Measurements made 
on additional |110] and [100] samples having dif- 








(fe in to Gauss” 


FIG. 1. Oscillatory component of the transverse 
magnetoresistance for current along [110]; curve A: 
Bop; at 1.7°K; curve B: Bog, at 4.2°K; curve C: Byijp 
at 1.7°K. pg, is the zero-field resistivity. 
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ferent carrier concentrations yielded oscillatory 
effects of lower amplitudes but leading to results 
in agreement with the data presented. Oscillatory 
effects in the longitudinal direction were not re- 
solved in these samples, the amplitudes being at 
least an order of magnitude smaller than the 
transverse effects. In all cases, magnetoresist- 
ance and Hall voltage measurements were made 
using standard four-probe techniques. Data were 
taken for all combinations of current and mag- 
netic field directions on two sides of the crystal 
and the results averaged. The Hall voltage and 
the magnetoresistance above 1 kgauss were 
linear functions of the magnetic field with a 
small oscillatory component superimposed. The 
oscillatory part was obtained by biasing out the 
linear term with a signal from a rotating-coil 
gaussmeter while slowly sweeping the magnetic 
field. 

Periodic oscillations in the magnetoresistance 
will occur as a result of the Landau levels pass- 
ing through the Fermi level EF. In a degenerate 


system, which strictly applies in our case, 


. (n h y? 7 


*, 2m 1 


0D 
where m))* =m7"*N,**K" is the reduced density- 
of-states effective mass, N, is the number of 
valleys, andK=m,*/m*, the mass anisotropy 
for a given valley. With E;/hw=a/B, where 

w =eB/m,* is the cyclotron frequency, the pe- 
riod, in units of 1/B, can then be expressed as 
1/a. The value of a is 


31.4 (n/N )**F(K), 


for B in units of 10* gauss and 7 in units of 
10'*/cm!%. 

For a (111) model, oscillations of a single pe- 
riod would be expected for B along a [100] direc- 
tion because of the single effective mass asso- 
ciated with this direction. For this orientation 
the experimentally determined value of a is 
4.9+ 0.3 as obtained in Fig. 2(a). The effective- 
mass ratio K determined from weak-field mag- 
netoresistance is 4.5+ 0.5 at 300°K and decreases 
somewhat at lower temperatures, presumably 
due to an admixture of ionized impurity scatter- 
ing. Using this value for K, and the expression 
for f(K),o, as given in Table I, eight- and four- 
valley models yield a =3.3+0.1 and 5.2+0.1, 
respectively. Thus a four-valley ellipsoid model 
is established. 

The treatment of transverse magnetoresistance 
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FIG. 2. (a) Nodes of the [001] and [110] oscillations 
in Fig. 1 plotted at half-integer points for the deter- 
mination of the periods and the phase shifts. (b) Decay 
of the Bop, oscillations vs 1/B at 1.7 and 4.2°K for the 
determination of the damping factors and thus the 
cyclotron effective mass (m,*)199. 


by Adams and Holstein‘ gives the following ex- 
pression for the oscillatory component for the 
case of at least several quantum levels and in- 
cluding only thermal broadening: 


(- a cos(2aME »/hw - $n) 


v i « 
a aw. a ED y 
p hwE 7] M=1 M* sinh(27°MkT/hw) ’ 


where p is the resistivity in the classical limit 





Table I. Values of m,* and f(K) for two B direc- 
tions and ellipsoids along the (111) axes. 





B f(K) 





m,.*KY*13/(K +2)” = K[3/(K +2)]”? 
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m, K K 


m,*K¥13/(2K +1)? — K"[3/(2K + 1))”? 
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and Pog is the oscillatory component. Our ex- 
perimental results indicate that the damping is 

so severe that only the first term in the summa- 
tion need be considered. This striking attenua- 
tion in the amplitudes for decreasing B can be 
attributed to both thermal broadening and broaden- 
ing of the Fermi level due to inhomogeneities in 
the crystal, the latter being the dominant effect. 
Collision broadening should be very small be- 
cause of the high mobility. 

For B in the [110] direction, two periods would 
normally be observed. However, because of the 
damping terms, the component due to the high 
effective mass will be very small after several 
periods, leaving predominantly the low-effective- 
mass oscillations. From the data, plotted in 
Fig. 2(a), the ratio a,,/a@,9, determines a value 
of K in approximate agreement with the value de- 
termined from weak-field magnetoresistance at 
300°K. In the [100] case, the phase agrees with 
that predicted by Adams and Holstein. The phase 
in the [110] case is somewhat shifted from the 
predicted value as might be expected from a 
small admixture of the high-mass period. 

The damping term which includes the effect of 
inhomogeneity and thermal broadening can be 
written as exp|-(8,+8,)/B], where 8, = 200(m,*/ 


m p*)(An/n) and g, = 14.8 m,*T. The observed 
damping, 8,>7, indicates a value for An/n of 


about 10% which is consistent with the measured 
inhomogeneity. 

To obtain the value of the total damping term 
[where the effect of inhomogeneity broadening is 
included in Eq. (1)], ApB~“*/p,T vs 1/B was plot- 
ted in Fig. 2(b). p in Eq. (1) was taken to be the 
actual resistivity at a given field, i.e., pxp,B; 
however, the actual magnetic field dependence 
in front of the oscillatory term has only a small 
effect on the determination of the damping factor. 
Since the difference in the slope at the two tem- 
peratures arises from the change in the value of 
the g, term, an effective mass m,* =0.035+ 0.005 
could be determined. Thus for K =4.5, m7* 
= 0.025+ 0.005, m* =0.12+0.03, and mp*= 0.10 
+0.02. Thermoelectric measurements at T < 10°K 
have given m))* = 0.08 + 0.03.? 

More complete results will be published in the 
future, including studies of longitudinal magneto- 
resistance and Hall oscillations. 

We wish to thank V. J. King of this laboratory 
for the preparation of the samples. 
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TRANSVERSE-EVEN VOLTAGE: A HIGH-FIELD GALVANOMAGNETIC EFFECT 
ASSOCIATED WITH OPEN ORBITS IN METALS 


J. R. Klauder and J. E. Kunzler 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received January 30, 1961) 


The magnetoresistance, determined from the 
voltage measured parallel to the current, repre- 
sents only one component of the even electric 
field, the even part of the electric field being 
that part which is unchanged by a reversal of the 
magnetic field. The other electric field compo- 
nents lie transverse to the current, and account, 
for example, for the well-known planar Hall 
effect. Here, we report observations in copper 
of such transverse voltages, and demonstrate 
a significant application of the transverse-even 
voltage in the determination of open-orbit direc- 
tions. In addition to providing information about 
the shape of the Fermi surface, we find that 
these observations vividly support the concept 


of open orbits.’ 

When the magnetic field lies in the plane trans- 
verse to the current the planar Hall effect van- 
ishes. However, an even voltage usually per- 
sists due to the presence of a slight longitudinal 
misalignment of the Hall contacts to the sample. 
Such voltages are of one sign because they repre- 
sent reproductions of the magnetoresistance on 
a reduced scale, and indeed, they are not true 
transverse voltages. In contrast, the transverse- 
even voltages which we have observed in high- 
purity copper single crystals for high magnetic 
fields transverse to the current do not scale with 
the magnetoresistance and characteristically 
exhibit changes of sign.” This new behavior, im- 
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plicit in the definitive work of Lifshitz and Pes- 
chanskii,® arises from the existence of a bulk 
anisotropic magnetoresistance tensor whose 
principal axes are sensitive to magnetic field 
orientation because of the presence of open orbits. 
These new voltages make it possible to determine 
open-orbit directions with significantly less ex- 
perimental data than would otherwise be required. 
In addition, these new voltages are more sensi- 
tive to some higher order* groups of open orbits 
than are magnetoresistance data alone. The ex- 
pected theoretical behavior of these new voltages 
is discussed under certain simplifying assump- 
tions. Predictions based on these arguments are 
compared with observed results. 

Consider those field directions at high mag- 
netic field strength such that an open Fermi 
surface leads to a non-negligible number of open 
orbits whose net direction is along a single direc- 
tion, say for convenience the x direction. In 
this case, there is only one dominant even term 
in the magnetoresistance tensor, p,, =H*a. The 
positive factor a is, to a first approximation, 
field insensitive, depending rather on the details 
of the scattering mechanism in a complicated 
manner. Let & be a unit vector in the direction 
of the open orbits andj be the current density; 
then 

E, =H7a2(8-)) (1) 
expresses the general relationship between the 
open-orbit direction, the current, and E,; the 
component of the electric field due to open orbits. 
The longitudinal voltage V;= D; -E, is proportional 
to the resistivity, while the two transverse-even 
voltages are defined by Vt, =D; Eo and V;, 
=D, Eo. In these expressions, the dD; are three 
orthogonal vectors, either parallel (i=/) or trans- 
verse (i =t,,¢,) to the current, whose lengths are 
proportional to their respective magnetoresis- 
tance- or Hall-probe separations. If, as the 
orientation of the magnetic field is changed, a 
new set of open orbits is encountered whose net 
direction is parallel to 2°, then a new set of 
voltages V7’, V;’, and V; ’ become significant. 
Various open-orbit directions may be encount- 
ered for a continuous rotation of the magnetic 
field. For each encounter the sign of V7 will be 
the same since V,<(D;-8). On the other hand, 
the uniformity of sign need not be maintained 
for the transverse voltage since V;; « (Dy, -%) 

x (Q- -D)), whose sign depends on the relative 
orientation of Q; this behavior accounts for the 
changes of sign of the transverse-even voltage 
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that are observed. As a consistency check on 
the presence of open orbits whose direction has 
been ascertained from previous measurements, 
it suffices to compare the predicted transverse- 
even voltages V; PF and V; PF, where 


Prey aya 
/* ed as” 2)/(D, $2), (2) 


with the corresponding observed values. 

The direction & of the open orbits, whose 
presence is signified by a relatively large value 
of Vj, is determined by 


Q = (normalization)[(V)/D)), (V, /D,)s (Vv, Pb 


(3) 
as expressed in the “sample” coordinate system. 
This procedure involves considerably less ex- 
perimental effort than the previous way in which 


Q was determined from magnetoresistance data 


alone, namely, by tracing out a sequence of 
points of high magnetoresistance in the plane 
transverse to 

One of the strong points of high-field magneto- 
resistance data is the information they can give 
regarding suitable topological models of the 
Fermi surface through the presence of variously 
directed sets of open orbits which the Fermi sur- 
face must support. For an understanding of the 
complex metals a qualitative model.of the Fermi 
surface should precede the “pinpointing” of cer- 
tain surface details. The use of the transverse- 
even voltages should more readily enable the 
Fermi surface to be determined. We now briefly 
discuss some of our transverse-even voltage ob- 
servations in copper. 

Figure 1(a) shows the magnetoresistance ob- 
served in copper at 18 kgauss in a transverse 
rotation about an axis near to the [511] direction. 
The peaks of magnetoresistance are proportional 
to H? and occur at orientations of the magnetic 
field for which a sizable number of unidirectional 
open orbits exist. In Fig. 2 we represent, ina 
stereographic projection, the data of Fig. 1(a) 
idealized to either H* behavior, represented by 
short heavy bars, or to saturation. In addition 
we show the three inequivalent types of principal 
magnetoresistance “ridges” —great circles on 
the unit sphere that have been found by earlier 
measurements to exhibit an H? behavior.*»* The 
interpretation of each magnetoresistance ridge 
requires the existence of open orbits directed 
normal to the plane of that ridge; thus, open 
orbits are directed along (100), (110), and (111) 
axes. Based on this knowledge of the open-orbit 





VoLuME 6, NUMBER 4 


PHYSICAL REVIEW LETTERS 


FEBRUARY 15, 1961 








32 ' 
COPPER 147-J-1 
R273° K 
R4.2° K 
T= 4.2°K 
H = 18 KGAUSS 
’ ‘ 


! 


~ 8000 























| 
| 















































VOLTAGE IN xv 


80 40 ° -40 -80 
MAGNETIC FIELD ANGLE, 0 


FIG. 1. The longitudinal voltage and the two trans- 
verse voltages observed vs angle in a transverse rota- 
tion when the current direction is near to a [511] axis. 
(a) The longitudinal voltage proportional to the mag- 
netoresistance. (b) and (c) Solid curve: the even part 
(i.e., less the Hall voltage) of the transverse voltages 
observed across contacts 2-2’ and 3-3’ as shown in 
Fig. 2. Dashed curve: that portion, similar to 
Fig. 1(a), of the solid curve which is due to contact 
misalignment. The difference of the solid and dashed 
curves represents the transverse-even voltage present 
because of the anisotropy of the magnetoresistance 
tensor due to open orbits. 
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FIG. 2. Stereographic representation of the mag- 
netoresistance data of Fig. 1(a) along with the known 
contours of H? behavior. The sequence of short heavy 
bars schematically represents the peaks of magneto- 
resistance observed in the transverse rotation. The 
great circular lines represent the three inequivalent 
types of magnetoresistance ‘ridges’ previously ob- 
served which indicate the presence of open orbits, 
along (100), (110), and (111) axes, perpendicular to 
the ridges of H* behavior. 


directions for copper, we have calculated the 
voltage in the two transverse directions that 
should be observed for each direction of the mag- 
netic field in the transverse plane such that open 
orbits exist. These results appear in column 3 
of Table I. 

The solid curves of Figs. 1(b) and 1(c) repre- 
sent the even part of the observed transverse 
voltages, while the dashed curves represent that 
component due to the misalignment voltage, 
which is a small fixed fraction of Fig. 1(a). In 
principle, the proportionality constant that re- 
lates each dashed curve to Fig. 1(a) can be de- 
termined from measurements at zero magnetic 
field; however, it turned out to be more suitable 
to determine these factors by a “best fit” to the 
data. Note the variation in sign; this fact is 
what first attracted our attention to these volt- 
ages. In column 4 of Table I we list the observed 
transverse voltages at various field angles which, 
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Table I. 
rotation. 


Predicted and observed values for the transverse-even voltages at various field angles in a transverse 





vy. pre 


ti 
Open-orbit 
direction 


Predicted® 
voltages (pv) 


obs 
Vv. : 


Open-orbit 
deviation 


angle 


i 
Observed* 
voltages (v) 





(110) 


(100) 


(110) 


(111) 
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(111) 


(110) 


. 6S 
- © 
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' 
= © 
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= 


47° (100), (111) 


0. 
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- 26 
- 49 
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- 99 


- 35 
-10 
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- 82 
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97 
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77 
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-23 
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-73 
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77 
22 


0.03 
-0.93 


-0.63 
1.09 


0.5 
-0.33 


of o of orf 
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-“ © eS 





@The top entry in each case is related to Fig. 1(b), while the bottom entry is related to Fig. 1(c). 
bThe top entry of 0.5 pv in this case was assumed to be correct in order to fix a composite open-orbit direction 


between [100] and [111]. 


according to Fig. 1(a), indicate the presence of 
open orbits so that Eq. (1) has validity. 

In column 5 of Table I we list the angle y, 
where cosy=Pre.Qtrue, the deviation between 
the true open-orbit direction and the direction 
predicted according to Eq. (3). In nearly every 
case, the angular deviation y does not provide 
an opportunity for confusion of the open-orbit 
direction with any other principal crystallograph- 
ic axis than the correct one. The simple analy- 
sis of Eq. (1) ignores additional even terms of 
lower order in the field; these become signifi- 
cant when 2-j =0. The vanishing of 2-j accounts 
for the absence of any peaks in Fig. 1 when 
@= -34°, and its smallness affects the accuracy 
of the predicted voltages at 6 = -78° and 6 =10°. 
The transverse voltages appear more sensitive 
[see field angles between -30° and -40° in Figs. 
1(b) and 1(c)] to the “little ridges” which have re- 
cently been observed in copper,*** and which may 
be expected in other metals, whose origins lie in 
the existence of small bands of higher order open 
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This direction was then utilized to predict the value -0.31 for the bottom entry by Eq. (2). 


orbits supported by the Fermi surface. 
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Roy. Soc. (London) A250, 325 (1957)], cyclotron reso- 
nance [e.g., D. N. Langenberg and T. W. Moore, Phys. 
Rev. Letters 3, 328 (1959)], and magnetoacoustic meas- 
urements [e.g., J. D. Gavenda and R. W. Morse, Bull. 
Am. Phys. Soc. 4, 463 (1959)] give less direct informa- 
tion concerning open orbits. 

*Similar voltages have been reported by J. Yahia and 
J. A. Marcus, Phys. Rev. 113, 137 (1959), for gallium, 
but are not interpreted. 

5], Lifshitz and V. Peschanskii, J. Exptl. Theoret. 
Phys. (U.S.S.R.) 35, 1251 (1958) [translation: Soviet 
Phys. —JETP 8, 875 (1959)). 
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53, R. Klauder and J. E. Kunzler, Proceedings of 
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ATTENUATION OF LONGITUDINAL ULTRASONIC VIBRATIONS BY 
SPIN- PHONON COUPLING IN RUBY 


E. B. Tucker 
General Electric Research Laboratory, Schenectady, New York 
(Received January 30, 1961) 


The attenuation of 9.1-kMc/sec longitudinal 
ultrasonic vibrations due to spin- phonon interac- 
tion has been observed in pink ruby. The value 
of the magnetoelastic coupling constant has been 
obtained from the experimental value of the at- 
tenuation. 

Calculations of the “direct process” spin lat- 
tice interaction predict that the coupling involved 
should be quadrupolar in character’~$ (i.e., 
through combinations of operators such as S,’, 
S,S, +S,S,, etc.) for non-S-state ions of the 
iron group such as Cr***. For longitudinal 
waves along the z axis, in this case both the 
trigonal axis of the ruby and the axis of quantiza- 
tion, the interaction may be represented by a 
perturbation of the form 

H’=Ge_S%, 
zzz 
where G is the magnetoelastic coupling constant 
and e,, is the elastic strain amplitude. The re- 
sulting transition probability for a spin to absorb 
a quantum of energy from the phonon field is 


Gre atv) \(1S,71) 1?/2h?. 


g(v) is the shape factor for the spin-lattice in- 
teraction, the maximum value having been pre- 
viously determined as ~10-* sec.‘ (|S,?|) repre- 
sents the matrix element of the spin operators 
between the levels involved. Experimental data 
for the angular variation of the transition pro- 
bability, determined using saturation techniques 
similar to those previously described,*~* are 
given in Fig. 1. Theoretical transition probabili- 
ties assuming quadrupole coupling have been 
computed using the ruby calculations of Chang 
and Siegman.” The experimental data agree 
quite well with the theoretical curves superposed 
on a background presumed to be due to trans- 
verse components, present by virtue of some 
longitudinal- transverse mode conversion in the 
crystal. At 6=0° the + 1/2 (high field 1-3) tran- 
sition is, as expected, not affected by the pho- 
uons while the 1-2 and the other 1-3 transition 
(both pure Am =+1) interact only with transverse 
waves. It has previously been reported that 
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FIG. 1. Angular variation of the ultrasonic transi- 
tion probability for longitudinal acoustic vibrations 
along the c axis of pink ruby. The transitions are 
labeled according to the energy levels involved, with 
the levels being numbered from the highest (1) to the 
lowest (4). 


Cr*** in ruby obeys the quadrupole selection 
rules.® 

The value of the attenuation per centimeter of 
path is 


a= Gg(v) ICIS, *|) Phy /(2h*pv*), 


where n is the population difference per centi- 
meter between the levels involved, v is the ve- 
locity of longitudinal acoustic waves, and p is 
the density of ruby. Since we appear to be justi- 
fied in taking the theoretical values for {|S,?|), 
the determination of a value for G? requires an 
experimental evaluation of a. 

The attenuation observation was made using a 
pulse echo technique similar to that used to de- 
termine the attenuation in quartz.*»’° A double- 
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cavity arrangement,’ ultrasonic drive and spec- 
trometer cavities, allowed the amplitude of the 
absorption signal to be monitored while the echo 
experiment was in progress. The ultrasonic 
drive cavity was pulsed at a rate of 10 cps using 
half- microsecond pulses of 90-watt peak power. 
A superheterodyne detection system was used to 
monitor the power from the drive cavity through 
a 3-db directional coupler. The crystals used 
were in the form of rods 3 mm in diameter, a 
13-cm length of X-cut quartz and a 1}-cm length 
of Linde pink ruby. The axis of the ruby rod is 
25 minutes from being parallel to the c axis. 
The ends of both quartz and ruby are optically 
polished flat, and are parallel to within one min- 
ute. The bonding material is indium. 

Figure 2 illustrates the echo pattern observed 
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FIG. 2. Echo patterns observed: upper trace with 
the magnetic field adjusted for resonance at @=60° on 
the 2-3 transition, and lower trace with the magnetic 
field just off the resonance value. The origins of the 
echoes are indicated between the traces. 
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at 1.5°K. The third echo has traversed the ruby 
rod four times for a total path length of 5 cm in 
ruby and 34 cm in quartz. A longer path length 
in the ruby would be desirable but, for the par- 
ticular rod lengths used, echo number four and 
all succeeding echoes are composites of reflec- 
tions from the bond and from the end of the ruby. 
No useful conclusions may be drawn from them. 
The difference in amplitude of the third echo in 
the top and bottom trace represents the absorp- 
tion of phonons by the spin system. The ampli- 
tude of the third echo averages 20% lower with 
the magnetic field at the value required for re- 
sonance at @=60° on the 2-3 line compared to its 
amplitude with the magnetic field just off the re- 
sonance value. This transition represents the 
optimum combination of experimental ultrasonic 
transition probability and Boltzmann factor. The 
conclusions below are based on a number of ob- 
servations made on this transition. Attenuations 
have also been observed for the 1-2 transition at 
@=20° and for the 3-4 transition at 6 = 70°. 

The attenuation observed is 0.2 for the 5-cm 
path or 0.04 per cm. The corresponding value 
of G? is 1.5x10°** erg’. If a direct process is 
assumed for the relaxation process and this 
value of the coupling constant is assumed to ap- 
ply to the transverse as well as to the longitudi- 
nal waves, the corresponding spin-lattice relaxa- 
tion time is 0.5 second, assuming Orbach’s for- 
mula" 


1/T, = G*1{|S,*1) wT /(2ah*pv'), 


with T =4.2°K and a velocity (the mean of the 
longitudinal and transverse velocities) of 9 x 10° 
cm/sec. There are experimental relaxation- 
time measurements™>** which agree quite well 
with this figure, while others**,’* are an order 
of magnitude lower. In any case the experimental 
relaxation time results do not show the variation 
with frequency predicted by this formula. 
Attenuations were not observable at @=0° for 
any of the three absorption lines allowed by the 
microwave selection rules. Theoretically this 
is what one expects since none of the transitions 
should be directly affected by a longitudinal wave. 
The sensitivity in the present experiment is, 
however, such that it is not reasonable to draw 
any conclusions concerning the relative propor- 
tions of longitudinal and transverse waves from 
these results for 6=0°. The value of the attenua- 
tion observed for coherent longitudinal waves 
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should be accurate to better than a factor of two. 
More accuracy should be obtainable with a judi- 
cious choice of rod lengths. 

It is a pleasure to acknowledge numerous dis- 
cussions with Dr. N. S. Shiren and Dr. E. H. 
Jacobsen, and the loan of the pulser from the 
latter. T. G. Kazyaka helped take a large part 
of the data. 
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ANOMALIES IN THE MAGNETIC SUSCEPTIBILITY AND ELECTRICAL RESISTIVITY OF VANADIUM 


Jean-Paul Burger 
Laboratoire Pierre-Weiss, Institut de Physique, Strasbourg, France 


and 


Michael Alan Taylor 
Solid State Physics Division, Atomic Energy Research Establishment, Harwell, England 
(Received January 30, 1961) 


An anomaly in the temperature dependence of 
the electrical resistivity of vanadium has been 
reported by several workers.’ The tempera- 
ture at which this anomaly occurs and its exact_ 
form seem to be markedly dependent upon 
specimen purity. The anomaly has been reported 


in the same specimen. Resistivity data on some 
vanadium-chromium and vanadium -cobalt alloys 
are also given. 

The susceptibility of vanadium between 20 and 
293°K has been discussed by Childs, Gardner, 
and Penfold.® Detailed measurements of the 









al as occurring at the following temperatures: variation of susceptibility with temperature have 

’ 200°K,* 210°K,? 190°K,° and 245°K.* Loomis and not been previously reported, although it is 
Carlson’ also report anomalies in the following known*’® that the susceptibility increases slight- 
physical properties at similar temperatures: ly with decreasing temperature. The variation 
lattice constant, Young’s modulus, internal fric- of susceptibility with temperature for a speci- 
tion, thermal expansion, and intensity of x-ray men of 99.9% purity Van Arkel vanadium kindly 

reflections. Hren and Wayman,” however, do not supplied by the Batelle Memorial Institute, 

e. find any discernible deviation from linearity in Columbus, Ohio is given in Fig. 1. The impur- 
a plot of Young’s modulus against temperature. ities (in wt. %) in the vanadium are as follows: 
The mechanism responsible for this behavior is C, 0.005; Si, 0.001; Cr, 0.001; Fe, 0.04; Al, 
not recognized. 0.005; Cu, 0.001; H, 0.001; N, 0.008; O, 0.020. 

The purpose of this communication is to report The measurements were taken on a Fo&x-Forrer 












translation balance in Strasbourg. The relative 
error associated with each susceptibility meas- 
urement is indicated on the graph. The absolute 


a previously undetected anomaly in the temper- 
ature dependence of the magnetic susceptibility 
and its association with the resistivity anomaly 
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FIG. 1. Variation of sus- 
ceptibility with temperature. 
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accuracy of the susceptibility measurements is 
+2%; the temperature measurements are accur- 
ate to+2°K. The results of three different runs 
are indicated (---, -+--, —). The susceptibility 
was independent of magnetic field. 

The variation of resistivity with temperature 
of the same specimen of vanadium is plotted in 
Fig. 2. The resistivity measurements were 
taken at Harwell using standard potentiometric 
techniques. The absolute value of the resistivity 


21+ 
3 
re) 
~ 
= 
= 
e) 
| 
© 20F 
a 
Y 
= 
VY 
> 
= 
> 19e 
— 
Ww) 
7) 
W 
« 

1 1 








i Nl i 
225 230 235 240 245 
TEMPERATURE (°K) 


FIG. 2. Variation of resistivity with temperature. 
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is only accurate to + 5%, but the relative accur- 
acies of the measurements are very good, and 
the measurements are perfectly reproducible. 
The resistivity at room temperature is 23+1 
pohm/cm. 

The suggestion by Rostoker and Yamamato* 
that the resistance anomaly is due to a trans- 
formation to a body-centered tetragonal phase 
stable at low temperatures has been shown to be 
incorrect. No such phase has been detected in 
low-temperature x-ray investigations by Loomis 
and Carlson' and Hren and Wayman’ or the pres- 
ent authors. Furthermore no well-defined dis- 
continuity, such as would be expected for a first- 
order phase change, occurs in the resistivity- 
temperature plot. 

Shull and Wilkinson’ report that neutron dif- 
fraction shows no detectable antiferromagnetic 
reflection at the (100) position at 20°K and that 
any moment present must be less than 0.1 p;. 
this conclusion has also been reached from pre- 
liminary neutron diffraction work at Harwell. 
However, the addition of 2% Cr into solid solution 
in vanadium lowers the temperature of the re- 
sistivity anomaly by 15°C, and the addition of 
5% Co by approximately the same amount; this 
is strongly reminiscent of the results of de Vreis* 
on dilute solutions of transition-metal elements 
in chromium and seems to indicate that the 
anomaly is, after all, magnetic in origin and is 
possibly due to the onset of antiferromagnetism. 
The fact that no antiferromagnetic moment has 
yet been observed in vanadium does not prove 
conclusively that antiferromagnetism does not 
exist. The time taken by a thermal neutron in 
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passing from one atom to the next is ~10~** sec, 
whereas the time taken by an electron near the 
Fermi surface is ~10~** sec; therefore if the 
antiferromagnetism is subjected to some process 
fluctuating more rapidly than 10~** sec, it is very 
probable that no (100) peak would be observed. 

The problem of the nature of the anomaly is 
being further investigated. 

Thanks are due to Dr. W. E. Gardner (Harwell) 
for valuable discussions and for providing the 


alloy specimens, and to Mr. T. E. Baker (Harwell) 


for carrying out the x-ray investigations. 
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CAPTURE OF SLOW NEUTRONS BY NUCLEI BOUND IN CRYSTALS* 


H. E. Jackson, L. M. Bollinger, and R. E. Coté 
Argonne National Laboratory, Argonne, Illinois 
(Received December 29, 1960) 


The intense interest in the recent discovery" 
of recoilless emission and absorption of y rays 
came as a sharp reminder of the classic paper 
in which Lamb foresaw the comparable effect for 
neutrons.” For neutrons the effect is manifested 
by the marked influence that crystalline binding 
of the target nucleus can have on the shape of a 
resonance. We have undertaken to determine 
whether the expected distortion in shape can be 
detected and, if so, to inquire whether it is large 
enough to give useful information about crystal 
dynamics. 

As was shown by Lamb, it is convenient to 
consider two sets of conditions under which the 
shape of a resonance observed with a crystalline 
target may be expected to depart from the shape 
that applies to a gas target at the same temper- 
ature. First, if the temperature T is not much 
greater than the Debye temperature @ and if k@ 
is much less than the classical recoil energy R, 
the shape remains that of a Doppler-broadened 
resonance for a gas target but the effective tem- 
perature needed to characterize the magnitude 
of the Doppler broadening is greater than the 
physical temperature of the sample. The ratio 
of the effective temperature 7’ to the true tem- 
perature is given by the equation® 


=. C,,(9/T) 36/T ' 6 


T 24 * 4[expo/T)-1) *7’ 


where C,,(@/T) is the Debye specific heat. Sec- 
ond, under the more extreme conditions [SR <k0é 





and T <6, the shape departs significantly from 
that of a resonance in a gas target at any tem- 
perature. The conditions R<k@ and T<6 are 
necessary for a significant fraction of the neu- 
tron captures to take place without recoil of the 
individual target nucleus, and the condition '<R 
is necessary so that the recoilless peak and other 
structure in the phonon-excitation spectrum of 
the crystal is not entirely obscured by the natural 
width of the nuclear level. 

In nature there are few materials that satisfy 
the conditions needed to give a significant de- 
parture from the elementary Doppler -broadened 
shape. Probably the most favorable case is the 
6.7-ev resonance in U***. For it T =0.026 ev and 
R =0.028 ev so that, if 6 is sufficiently high, the 
predicted effects should be observable. Con- 
sequently we have studied the line shapes for 
capture in two uranium samples, one a metallic 
foil of normal uranium, for which ké@=0.017 ev, 
and the other a sample of U,O,, for which ké 
= 0.038 ev at liquid nitrogen temperature, accord- 
ing to the available data on specific heats. 

The shape of the uranium resonance was de- 
termined by transmission measurements with 
the Argonne fast chopper* operated under condi- 
tions such that the time-of-flight resolution width 
was about 0.038 usec/m (full width at half max- 
imum), a value which corresponds to an energy 
spread of 0.018 ev at 6.7 ev. Both the metallic 
foil and the oxide sample were studied at room 
temperature and at liquid nitrogen temperature, 
and the foil was studied at the temperature of 
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liquid helium. For each of these conditions of 
operation, the neutron transmissions of the 
sample were measured over the energy range 
from about 4.5 to 30 ev in a series of several 
independent runs. The background counting rate 
was continuously monitored by passing the neu- 
tron beam through filters of platinum and silver 
that were thick enough to remove all timed neu- 
trons in the neighborhood of the resonances at 
11.9 and 5.2 ev. 

Some of the data obtained are presented in 
Fig. 1. The resonance shape for the metallic 
foil at room temperature was compared with the 
gas-model shape of a Doppler-broadened reso- 
nance by curve fitting. Since crystalline binding 
is not expected to play a significant role for 
these data, it was not surprising to find a good 
agreement, as indicated by a x” test, between 
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FIG. 1. Total cross-section data for the neutron 
resonance at 6.65 ev in uranium metal at temperatures 
of 297°K and 4.3°K. The solid curves are the shapes 
predicted by the gas model with the physical sample 
temperature in (a) and with the effective sample tem- 
perature (75°K) in (b). The dashed curve of (b) is the 
shape predicted by a Debye crystal model in which 
k@=0.017 ev. 
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the experimental and calculated shapes. The 
small systematic deviation of the data from the 
curve in the region of high cross section is due 

in part to the finite resolution of the chopper, 
which was not included in calculations of the line 
shapes. The parameters obtained by curve fitting 
and area analysis of the data are energy=E, 
=6.652 ev, peak cross section =o, =20700+ 1200 
barns, I =0.0265+ 0.0015 ev, and o,I =547+10 

ev barns. These values are in satisfactory agree- 
ment with previously published® results. 

Having determined reliable values of the reso- 
nance parameters from the data for the warm 
metallic foil, departures of the other data from 
the simple Doppler-broadened shape may now be 
detected by comparing the data with the shape 
calculated from the gas model at the appropriate 
temperature. Measurements of the metal at 
T =97°K are in fair agreement with the gas model. 
However, in contrast to this, the data on the 
metal at T =4.3°K (Fig. 1) and those on the oxide 
(Fig. 2) at both T =297°K and 97°K are not even 
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FIG. 2. Total cross-section data for the resonance 
at 6.65 ev in U,;O, at temperatures of 297°K and 97°K. 
The solid curves are the shapes predicted by the gas 
model with the effective sample temperatures [334°K 
in (a) and 181°K in (b)]. The dashed curve of (b) is the 
shape predicted by a Debye crystal model in which 
k@=0.038 ev. 
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qualitatively consistent with the gas-model shapes 
corresponding to the physical temperature. Thus 
we proceed to compare the data with the gas- 
model shapes corresponding to the effective 
temperatures of the samples. This kind of a 
comparison is made in Figs. 1(b), 2(a), and 2(b), 
where the solid curves are the Doppler-broadened 
resonance shapes for effective temperatures 75°, 
334°, and 181°K, respectively. The data for the 
metallic foil at liquid helium temperatures agree 
with the calculated curve, insofar as the maxi- 
mum cross section and general shape of the reso- 
nance are concerned. However, the energy of 

the maximum is observed to be displaced toward 
lower energies. This displacement cannot be an 
instrumental effect resulting from an instability 
in the zero time of the neutron time-of-flight 
system since, on comparing the data for the 
warm and cold metal foil, the position of the 6.7- 
ev resonance is observed to move with respect 

to the other resonances in U*** and the resonances 
in platinum and silver. Thus the displacement 

of the maximum is believed to be real. It pre- 
sumably is an indication that an abnormally 

large fraction of the neutrons are captured with 
small energy transfers to the crystal lattice. 

For these events, the incident energy required 

to excite the nuclear level is less than that re- 
quired for free nuclei. 

Because of its high Debye temperature, we 
would expect to observe much larger departures 
from the gas-model shape for a target of U,O,. 
This expectation is fulfilled, as is shown by the 
data of Figs. 2(a) and 2(b). The line shape at 
room temperature is adequately described by 
using the effective-temperature approximation. 
However, at liquid nitrogen temperature the 








data show a marked departure from this approx- 
imation. 

To understand the departure of the observed 
resonance shape from the effective temperature 
approximation one must assume a specific model 
of phonon excitation. Assuming the Debye model 
to be applicable, we have calculated the shape 
using the asymptotic expansion of Nelkin and 
Parks.® For the metallic foil the calculated re- 
sult, as given by the dashed curve of Fig. 1(b), 
is seen to be in excellent agreement with the ex- 
perimental data. For the oxide sample, however, 
the calculated curve is inconsistent with the data 
of Fig. 2(b). This discrepancy is not surprising 
since the thermodynamic behavior of U,O, at low 
temperature indicates that the Debye phonon spec- 
trum may not be appropriate. An analysis is now 


in progress to attempt to interpret the observed 
shape of the resonance in terms of a more real- 


istic model. 
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ANTISHIELDING OF NUCLEAR ELECTRIC HEXADECAPOLE MOMENTS* 


R. M. Sternheimer 
Brookhaven National Laboratory, Upton, New York 
(Received January 19, 1961) 


Up to the present time, evidence for a nuclear 
electric hexadecapole (16-pole) moment (HDM) 
seems to have been obtained only in one experi- 
ment,’ probably on account of the fact that the 
effects of the HDM are masked by the much lar- 
ger effects of the nuclear quadrupole moment Q. 
Nevertheless, it is expected that nuclei with spin 
J 22 will have a nonvanishing HDM, and the 
measurement of its value would obviously con- 
stitute a check on any proposed nuclear wave 
function from which the HDM could be calculated. 
The purpose of this Letter is to point out that 
for the case of ions, the effect of the nuclear 
HDM will be considerably amplified by antishield- 
ing effects of the same type as have been calcu- 
lated*~* and observed® for nuclear quadrupole 
moments. An estimate will be made for the anti- 
shielding factor n.. for the HDM. Here the HDM 
which is induced in the closed electron shells is 
written as® Hj,q=-NeoH, where H is the nuclear 
HDM. The antishielding effect is due to the ex- 
citation of electrons from the closed d shells 
into higher d states, and similar nf +f perturba- 
tions, produced by the field of the nuclear H. 

The total HDM of the ion: 

Hon 7H +H ing = (1-0, JH, (1) 
interacts with the fourth-derivative terms of the 
potential of the lattice, V, evaluated at the nucleus; 
e.g., 8*V/az*, a*V/ax7az”, a*V/ay"az”. The aver- 
age value of the field of H is zero for s and p 
electrons, so that the antishielding effects will 
be present only for medium and heavy ions with 
closed d (and possibly f ) shells. 

We can define the nuclear HDM as follows: 


H=)) (r.*(35 cos*@. -30 cos"é , +3)) 
(p) ? 


where 7, is the distance of a proton from the 
center of the nucleus, and the angle 6, is meas- 
ured with respect to the z axis (direction of 
quantization). The expectation value in Eq. (2) 

is taken over the nuclear wave function with mag- 
netic quantum number m;=I, and the sum ex- 
tends over all protons in the nucleus. The poten- 
tial due to H acting on an electron is 


V.. =eHP,(8)/(8r*), (3) 


a 
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where ,¢ is the distance of the electron from the 
center of the nucleus, P, is the Legendre poly- 
nomial, and @ is the angle between the direction 
of the nuclear spin and the position vector of the 
electron. The perturbation of the potential ener- 
gy of the electron in rydberg units is given by 

U,, = -AP,(8)/(47°) ry, (4) 
where ¢ is in units of the Bohr radius ay, and H 
is in units az*. 

Upon proceeding in the same manner as in 
reference 2, one finds that the induced moment 
density due to the nd +d and nf +f (radial) modes 
of excitation is given by 


(r)dr =Hr* [$3 s2iMg "(adda yd 4) 


i, rad 


+ Lin, Op ut Ht “PI: (5) 


where the sums extend over the closed d and f 
shells of the core; u,’ is ry times the unperturbed 
(nd or nf) radial wave function, normalized ac- 
cording to [5 u,"dr=1; u’ 1,H 18 7 times the 
radial part of the perturbation, and is determined 
by the equations 


[ o Ave Bee a-2) 
-ustad| 3 Ce) |) © 
nl P| 
j- ugielut {l= Dar =0, (7) 


where V, and E, are the unperturbed potential 
and energy eigenvalues for the mi shell, and 
(1/r*) nj is the average value of 1/r° for the wave 
function u,’(ml). The coefficients 80/63 and 
112/99 in Eq. (5) result from the integration over 
the angular parts of the wave functions, and the 
summation over the magnetic substates m. These 
coefficients _ ~ 1) are obtained from 


C,(nl = l)=8 iz [fi1e"e, sinoas]”, (8) 
=-] 0 
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where 0/ is the spherical harmonic pertaining 
to magnetic quantum number m. From Eq. (5), 
the term 7.,.(nd +d) pertaining to the excitation of 
the nd shell is given by 


nid ~d)=(80/63){ “u,"ndu’, (nd ~ d)r*dr. 
, (9) 


An estimate of the order of magnitude of 
1.("l + l) can be obtained by using hydrogenic 
wave functions. This assumes that V, and E, 
can be written as V,=-2Z.2/r and E,=-Ze"/n® 
(rydberg units), where Z, is an appropriate 
effective Z. In this case, the problem of find- 
ing the perturbations can be solved exactly, in 
the same manner as in Eqs. (14)-(20) of refer- 
ence 3. This calculation has been carried out 
for 3d~d, 4d+d, and 4f~f, with the following 
results: 1.(3d +d) = -9274/(63Z 2) = -147.2/Z,; 
Noo(4d + d) = -105 892 /(63Z ,) = -1680.8/Z,; and 
Nol 4f +f ) = -24 524/(495Z 2) = -49.54/Z,. In order 
to evaluate these results, we must obtain the 
appropriate value of Z,. Z, will be a function 
of the nuclear charge Z; e.g., for 3d~+d, the 
value of Z, will depend on whether the 3d shell 
is the outermost shell of the ion, or whether 4s, 
4p, and 4d shells are also present, so that the 
3d electrons are essentially internal and see a 
considerable part of the nuclear potential -2Z /r. 

For an estimate of Z,, we may be guided to 
some extent by the corresponding results ob- 
tained for the quadrupole antishielding factor 
Yoo(nl+1). We have found, using hydrogenic 
wave functions®: y,.(2p + p)=-10.72/Z.; y..(3p = p) 
= -40.32/Z,; y..(3d ~d)=-7.77/Ze. We can equate 
these results for y,,.(m1 +1) with those obtained 
from Hartree-Fock wave-function calculations‘ 
to deduce the corresponding values of Z,. Thus 
for Nat, the wave-function result* y_(2p + p) 
=-5.16 leads to Z, =2.08; for Cl”, y,.(3p + p) 
=-56.5 gives Z,=0.71; for Cu’, y,,(3d +d) =-8.5 
gives Z,=0.91. It thus appears that if the mp or 
nd sheil considered is the outermost shell of the 
ion, the effective Z, is ~1-2. With Z,=1, one 
obtains n,,(3d +d) ~ -150; n..(4d +d) ~ -1700; and 
No(4f ~f)~-50. The values of ,,(nd~d) are 
very large, and may make it possible to observe 
the electric hexadecapole moment in quadrupole 
resonance spectra, using ions in crystals or in 
polar molecules. 

In addition to the strong antishielding effects 
due to the radial perturbations nd ~d and nf +f, 
one expects some shielding of the nuclear H due 
to the angular excitations of the core: ns ~g; 





np~f, np~h; nd~+g, andnd-+i. This shielding 
effect can be treated by means of the Thomas- 
Fermi model of the atom, in the same manner 
as was done previously’ for the angular part 
Qi,ang of the induced quadrupole moment. One 
obtains for the density of the induced HDM: 


a 12 
Fy ang *)4 = gH (xx)"dx, (10) 
where x and x are the Thomas-Fermi function 

and variable, respectively; x =(Z“*/0.8853)(r/ay). 
We note that Eq. (10) involves the same function 
(xx) as Qi,ang(), the only difference being 

that the coefficient? 3/10 is replaced by 1/6. [In 
general, the coefficient is given by (3/2)/(2L +1) 
for a nuclear 24-pole moment. Thus if the nu- 
cleus had an electric dipole moment D, the den- 
sity of induced dipole moment would be $D(xx)*dx. ] 
In view of Eq. (10), the total induced HDM is given 
by 


iy ( ~(yry2 
Hi ang" 2 f (xx)ax. (11) 
As was discussed in reference 2, this integral 
diverges logarithmically at large x, if the 
Thomas-Fermi function x (without exchange) is 
used. However, in the region of r<4ay, the 
Thomas-Fermi result for Qi, (r) agrees rea- 
sonably well with that obtained from wave-func- 
tion calculations, although (as expected) it 
smoothes over the maxima and minima of 
Qi,ang(”) which correspond to the contributions 
of the various electron shells. From wave-func- 
tion calculations, we have previously obtained 
the following values*’’ of yang: 0-60 for Na*, 
1.05 for K*, 1.4 for Cl”, and 2.9 for Cs*. If, 
following Eq. (11), one assumes that Neo, /Yeo,ang 
=(1/6)/(3/10)=5/9, one finds: Neo, ang = 9. 3 for 
Na*, 0.58 for K*, 0.78 for Cl, and 1.6 for Cs*. 
These results show that Neo, is of the order 
of 1, and is therefore very small compared to 
the contribution to n,. from the radial modes for 
atoms with closed d shells. Thus for medium 
and heavy atoms, 7, is determined predominant - 
ly by the contribution of the radial modes of ex- 
citation, which give rise to the antishielding. 

The preceding discussion pertains to the case 
of ions in crystals or in polar molecules, which 
appear to be the most hopeful possibilities for 
detecting the nuclear HDM. For the hyperfine 
structure of valence electrons in neutral atoms, 
the situation is somewhat different. In this case, 
the antishielding modes will probably not be 
nearly as important, since the corresponding per- 
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turbation of the electron density is located main- 
ly at rather large distances from the nucleus, 
whereas the interaction with the valence electron 
takes place very close to the nucleus on account 
of the 1/r° dependence of the potential Vj. The 
shielding modes make a relatively large contri- 
bution, since the induced density H; ang *) is 
located somewhat closer to the nucleus. In view 
of the results of reference 2, one may expect a 
shielding of the interaction of H with the valence 
electron of the order of ~10% for light atoms. 
For heavy atoms, with closed d shells, there 
may be a net antishielding due to the presence of 
the radial modes. 

Concerning possible experiments with ions in 
crystals, it seems that the following nuclei offer 
some good possibilities for observing the nuclear 
HDM: (1) Zn® ([=5/2, 3d*°4s"); Ge™ (1=9/2, 
3d'°4s74p"); (2) In™® (1=9/2, 4d°°5s75p); Sb*** 
(1=5/2, 4d*°5s*5p*); Sb’** (7=7/2); and (3) Bi?°® 
(1=9/2, 5d*°6s*6p*). For each nucleus, we have 
listed the nuclear spin J and the electron con- 
figuration of the outer shells. For all cases, 
the number of electrons outside a closed d shell 
is small (<5), so that the d shell is essentially 


external and Z, is small. The cases listed under 
(1), (2), and (3) correspond to a newly filled 3d, 
4d, and 5d shell, respectively. 

I wish to thank Professor E. M. Purcell for 
helpful discussions. 





*Work performed under the auspices of the U. S. 
Atomic Energy Commission. 

'T. C. Wang, Phys. Rev. 99, 566 (1955). 

*R. M. Sternheimer, Phys. Rev. 80, 102 (1950); 
84, 244 (1951); 95, 736 (1954); 105, 158 (1957). 

H. M. Foley, R. M. Sternheimer, and D. Tycko, 
Phys. Rev. 93, 734 (1954). 

‘R. M. Sternheimer and H. M. Foley, Phys. Rev. 
102, 731 (1956). 

See, for example, M. H. Cohen and F. Reif, in 
Solid State Physics, edited by F. Seitz and D. Turn- 
bull (Academic Press, New York, 1957), Vol. 5, 

p. 321; T. P. Das and E. L. Hahn, Nuclear Quad- 
rupole Resonance Spectroscopy (Academic Press, 
New York, 1958). 

®We use the same convention of signs as for the 
quadrupole antishielding factor y.; i.e., positive 
1. corresponds to shielding, while negative n,, cor- 
responds to antishielding of the nuclear HDM. 

"R. M. Sternheimer, Phys. Rev. 115, 1198 (1959). 














PION PRODUCTION BRANCHING RATIOS THROUGH THE 900-Mev RESONANCE* 


B. A. Munir,! E. Pickup,} D. K. Robinson, and E. O. Salant 
Brookhaven National Laboratory, Upton, New York 
(Received January 30, 1961) 


Carruthers and Bethe’s discussion! of the in- 
fluence of a 1-7 interaction on the 900- Mev pion- 
nucleon resonance involves the experimental 
branching ratio for the reactions 


1 +p—-n +n* +n, (1) 
and 
1 +p—1 +7°+p. (2) 


In the course of studying pion- proton interactions 
around 1 Bev, we have measured this ratio for 
an incident pion energy of 905+ 2 Mev. 

The reactions were produced by magnetically 
analyzed pions traversing the Brookhaven 20- 
inch hydrogen bubble chamber, in a field of 
12300 gauss. Events originating in a well-de- 
fined central portion of the chamber were meas- 
ured with digitized microscopes. The events 
were analyzed with a 704 computer, using the 
BNL “Tred” program for data reduction and the 
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“Guts” event-fitting program for kinematic iden- 
tification. The chamber was run so that it was 
possible to distinguish visually between protons 
and pions. Of 985 measured two-pronged events, 
530 were identified as elastic scattering,? 427 as 
single pion production, and 28 as double produc- 
tion events. 

Denoting by m, and n, the numbers of events of 
reactions 1 and 2, respectively, the branching 
ratio y was found to be 


vr =n,/n, = 269/158 = 1.702 0.17. 


Values of r from hydrogen bubble chamber ob- 
servations at energies E near 900 Mev are taken 
from various authors and are compiled in Table I. 
It is seen that the measurements are consistent 
with an approximately constant branching ratio in 
the 800-1100 Mev energy range. 

Early measurements’ of the cross sections of 
reactions (1) and (2) produced by 900- Mev pions 
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Table I. 





Branching ratios y for incident pion kinetic energies E. 









E (Mev) 810° 905 





960° 960° 1020° 


1100° 





Y 2.5+0.4 1.70+0.17 


1.25 +0.16 


1.56 + 0.19 2.0+0.5 1.67 +0.13 





4],, Baggett, Lawrence Radiation Laboratory Report UCRL-8302, 1958 (unpublished); L. Baggett and 
B. McCormack, 1958 Annual International Conference on High-Energy Physics at CERN, edited by 





B. Ferretti (CERN Scientific Information Service, Geneva, 1958), p. 68. 
by. Alles-Borelli, S. Bergia, E. P. Ferreira, and P. Waloschek, Nuovo cimento 14, 211 (1959). 
CE. Pickup, F. Ayer, and E. O. Salant, Phys. Rev. Letters 5, 161 (1960). 


4], Derado and N. Schmitz, Phys. Rev. 118, 309 (1960). 


©E. Pickup, D. K. Robinson, and E. O. Salant (to be published). 


in emulsions indicated r=0.7+ 0.3. However, the 
uncertainties involved may not rule out consist- 
ency with the bubble-chamber value. 

The Carruthers-Bethe model for the reactions 
(1) and (2), considering the amplitudes a,, a,, 
a, of final pion- pion isospin states ¢=0, 1, 2, 
predicts that 


la, |?/la,1? = (r - 3)", 


if the T=3/2 state is excluded (so that a,=0). By 
inserting the low emulsion value of r at 900 Mev, 
a sharp maximum in |a,|?/|a,|* at that energy 
was obtained, indicative of a resonant ¢=1 state. 
Clearly the new measurement of ry at 900 Mev 
contradicts this conclusion. However, it does 
not necessarily follow that some other model 
would not give a dominant ¢=1 state, even with 
the branching ratios observed. 


We take pleasure in acknowledging our indebt- 
edness to Dr. F. R. Eisler, Dr. N. I. Samios, 
and Dr. M. Schwartz for preparation of the pion 
beam, to the BNL Bubble Chamber Group, to 
the Cosmotron crew, and particularly to our 
scanners for their unstinting work. 





*Work performed under the auspices of the U. S. 
Atomic Energy Commission. 

Ton leave of absence from the Florida State Univer- 
sity, Tallahassee, Florida. 

ton leave of absence from the National Research 
Council, Ottawa, Ontario, Canada. 

‘Pp. Carruthers and H. A. Bethe, Phys. Rev. Let- 
ters 4, 536 (1960). 

This value is not corrected for scanning losses at 
small angles. 

3w. D. Walker, F. Hushfar, and W. D. Shephard, 
Phys. Rev. 104, 526 (1956). 





ELECTRONS IN THE PRIMARY COSMIC RADIATION* 


Peter Meyer and Rochus Vogt 
Enrico Fermi Institute for Nuclear Studies and Department of Physics, 
University of Chicago, Chicago, Illinois 
(Received January 16, 1961) 


The question of the presence of electrons in 
the primary cosmic radiation has frequently been 
discussed in the literature.’ Attempts to identify 
such a component have so far been negative. 
Critchfield, Ney, and Oleksa? quote an upper 
limit of 0.076 (min cm? sr)~! for electrons with 
energies exceeding 1 Bev. In this note we re- 
port an experiment which forces us to conclude 
that there exists a flux of primary electrons 





which enters the top of the atmosphere. 

Three balloon flights were carried out in Au- 
gust and September, 1960 from Fort Churchill, 
Manitoba which has a calculated geomagnetic 
cutoff rigidity for vertically incident particles of 
about 100 Mv (Quenby and Webber*). Data were 
obtained in each flight during the ascent and for 
approximately 10 hours under 3 to 5 g/cm? of 
residual atmosphere. A cross section through 
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64cm 






































FIG. 1. Cross section of the detector system. 


the detector is shown in Fig. 1. The telescope 
counter and the top Nal counter define the solid 
angle of acceptance for vertically incident par- 
ticles. Their range in lead is measured by de- 
termining the number of range counters that 
were triggered in sequence by at least one ion- 
izing particle. The pulse heights from the Nal 
counters I and II are measured and give the en- 
ergy loss of particles prior to entering the lead 
absorber and the energy loss of those particles 
able to penetrate the entire absorber thickness. 
This method discriminates against all upward 
moving charged particles except for those that 
have minimum energy loss in both Nal counters. 
The detector system is surrounded on four sides 


by plastic scintillation counters in anticoincidence. 


Protons which come to rest in the lead through 
ionization are readily identified by their range 
and the corresponding energy loss in Counter I. 
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In addition to these protons (and a particles) 
many events are observed which have minimum 
ionization loss in Counter I but end within the 
stack of lead absorbers. These events cannot be 
attributed to protons losing their energy by ion- 
ization only, since they would penetrate to count- 
er II, but may be due to (a) high-energy protons 
making a nuclear collision, (b) mesons, or (c) 
electrons entering the detector and producing a 
soft shower in the lead plates. 

The altitude dependence of those events is 
shown in Fig. 2. Near the transition maximum 
the intensity follows the altitude dependence for 
secondary electrons (dashed line), while at high 
altitude it begins to level off or even increase. 
The dashed line gives the approximate altitude 
dependence of the secondary electron intensity 
in air originating from 1 decay and was arrived 
at by assuming that the number of 7 mesons 
created at any depth in the upper atmosphere is 
proportional to the flux of high-energy protons 
present at this depth. It is normalized to the ob- 
served intensity at the transition maximum. Ob- 
viously, the observed altitude dependence cannot 
be explained on the basis of secondary electrons. 

In Fig. 3 we plot the observed differential range 
spectrum of the events under discussion at 3 to 4 
g/cm? atmospheric depth. The dashed line gives 























joo 
Atmospheric Depth (g/cm) 


FIG. 2. Intensity versus atmospheric depth. 
—— Mimimum- ionizing particles with range between 
0 and 120 g/cm? of lead (error limits shown are stand- 
ard deviations). - - - Secondary electrons (computed). 
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FIG. 3. Differential range spectrum of mini- 


mum-ionizing particles in lead under 3 to 4 g/cm? of 
air (error limits are standard deviations). - - - Pos- 
sible contribution by protons and mesons. 


the range spectrum that could be caused by pro- 
tons and pions under assumptions (a) and (b). A 
comparison of the curves shows that these par- 
ticles cannot alone produce the observed range 
spectrum. Therefore, a substantial fraction of 
these events is attributed to primary electrons. 
We then estimate the energy and flux of the 
electrons which appear to enter the top of the 
atmosphere. In particular, we wish to find out 
whether a portion of the electrons has energies 
exceeding the geomagnetic cutoff at Fort Church- 
ill for vertical incidence and must, therefore, 
be of interplanetary or galactic origin. Particles 
below the geomagnetic cutoff may originate as 
secondaries from the southern hemisphere. The 
electrons which enter the lead absorber will in 
general produce soft showers. We shall have to 
relate the “range” of the shower to the energy of 


the incident electron. Greisen and Thom‘ have, 
for several electron energies, determined the 
average range in lead at which all ionizing tracks 
of a shower have disappeared. From their data 
and Monte Carlo calculations of Wilson® we have 
derived the energy scale in Fig. 3. 

We see that primary electrons above the calcu- 
lated geomagnetic cutoff of 100 Mv are present. 
All data presented are from a flight on September 
8, 1960, at a time of normal solar activity. Meas- 
urements on other days lead to the same conclu- 
sions. 

We here restrict ourselves to deducing a low- 
er and upper limit for the primary electron flux. 
The lower limit must be derived by subtracting 
the largest possible contribution of secondary 
electrons, stopping protons, and mesons. This 
contribution was obtained by making the extreme 
assumption that in the range interval with lowest 
intensity, all the events are due to protons and 
mesons. The dashed line in Fig. 3, which gives 
the calculated proton and meson spectrum, is 
normalized at that point. The difference between 
the two curves, therefore, represents a lower 
limit for the electron flux. The possible contri- 
bution of secondary electrons to this flux was ob- 
tained from Fig. 2 and also subtracted. For rea- 
sons not discussed here the effect of upward 
moving neutrons and y rays is estimated to be 
less than 10% of the lower limits. 

The upper limits for the electron flux are given 
by the observed values without subtraction of a 
possible contribution by protons, mesons, or 
secondary electrons. This leads to the limits 
for the vertically incident electron flux that are 
shown in Table I. 

These are preliminary results and detailed 
calculations of the electron flux and energy spec- 
trum, its time variations, and the implications 
to the problem of the origin and modulation of 
cosmic radiation will be published elsewhere. 

For many contributions we are indebted to Dr. 
J. A. Simpson, Dr. E. N. Parker, T. Burdick, 

R. Weissman, and G. Lentz. 


Table I. Flux of vertically incident electrons (cm? sec sr)~', 





25 < E <100 (Mev) 


100 <E <1300 (Mev) E > 1300 (Mev) 





28 x 107% 
31x 1073 


Lower limit 
Upper limit 


3.5x107 0 
11x 107 8x 107 
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The series of balloon flights from Fort Church- 
ill was made possible through the cooperation of 
the National Research Council of Canada. 
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tary Particles and Cosmic Ray Physics (North Holland 
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337. 
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‘K. Greisen and H. Thom (unpublished). (We wish 
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able to us prior to publication.) 
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DIRECT = DECAY AND MUONIUM- ANTIMUONIUM TRANSITIONS 


Sheldon Lee Glashow 
California Institute of Technology, Pasadena, California 
(Received January 13, 1961) 


It now seems certain that the mass difference 
between K,,,.° and K,,,° is quite small (|5m |<10-° 
ev with 95% confidence’). Okun’ and Pontecorvo? 
discuss a connection between 6m and the exist- 
ence of AS =2 weak interactions. They show how 
such an interaction, if as strong as the AS=1 
weak interactions, might lead to a mass splitting 
of the order of electron volts. Indeed, if the 
AS =2 interaction were absent, a mass difference 
between the neutral K mesons arises only at 
second order in the weak interaction and the 
theoretical estimate, 5m =10-° ev, agrees with 
experiment. But suppose L,, (the weak-interac- 
tion Lagrangian) includes a AS =2 term (which 
would permit direct decays of cascade hyperons, 
i.ec., Z -n+n, =°+n+n°, and =°+p+7). The 
first order mass splitting, 5m, is equal to 
twice the matrix element, M=(K°|L,,|K°). Noting 
that M=(K°|L,, |K°) =(K°|C“*LC |K°), where 
C is the charge-conjugation operator; we find 
that only the part of the AS =2 coupling which is 
even under C contributes to 6m. Hence, a 
AS = 2 weak interaction which is odd under C does 
not lead to a large K,,.°, K,,° mass splitting. The 
small measured value of 6m implies nothing about 
the AS=2, C-odd interaction, but only requires 
that the AS=2, C-even interaction is no more 
than 10°° as strong as AS =1 weak interactions. 

Assume that a AS =2 interaction odd under C 
exists with similar strength as ordinary AS=1 
weak interactions [the four-Fermion coupling, 
G(Ay,”) (Ay, 75”) + H.c., is one example*]. From 
CP invariance it follows that any such interaction 
is also odd under P. Consequently, the effective 
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interaction responsible for direct decay of = into 
a nucleon and a pion is invariant under space re- 
flection for some assignment of relative =, n 
parity. The absence of a large 6m, although 
failing to forbid direct = decay, implies that such 
modes (if they occur at all) cannot display asym- 
metries (i.e., they “conserve parity”). No direct 
decays of = have yet been seen. However, exist- 
ing experimental evidence that direct decay is 
absent is exceedingly slight. 

An identical argument relates the appearance 
of muonium-antimuonium transitions to the possi- 
ble existence of interactions permitting e” +e-— 
u-+yu"." At first sight, looking for spontaneous 
transitions between muonium and antimuonium 
seems a more sensitive experiment than the di- 
rect search (utilizing, for example, the clashed 
electron beams soon to be available at Stanford 
University). But transitions between the 1S 
states of muonium and antimuonium are gene- 
rated only by that part of the (e° +e" ~ 1+ p per- 
mitting) interaction which is even under C. To 
prove this, merely read for ( K°| and { K°| in the 
first paragraph, the 1S state of muonium and of 
antimuonium. The interaction (fy,e)(iy,7<2) 
+H.c., odd under C, could hardly be detected by 
looking for the muonium-antimuonium transitions 
it fails to induce. It could be found directly (by 
charge-exchange scattering of muons on elec- 
trons, or bye +e °~+~yu~+y"~), or possibly, by 
searching for nondegenerate transitions between 
states of muonium and antimuonium in the pre- 
sence of external electromagnetic fields. 

We wish to thank Professor M. Gell- Mann for 
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his comments. ary phases (e.g., real K-meson parity) and L,, is odd 
under either C or P. 
'F, Muller et al., Phys. Rev. Letters 4, 418 (1960). ‘This coupling arises in a partially-symmetric model 
21. Okun’ and B. Pontecorvo, J. Exptl. Theoret. of weak interactions. S. Glashow (to be published). 
Phys. (U.S.S.R.) 32, 1587 (1957) [translation: Soviet 5B. Pontecorvo, J. Exptl. Theoret. Phys. (U.S.S.R.) 
Phys. -JETP 5, 1297 (1957)]. 33, 549 (1957) [translation: Soviet Phys. —JETP 6, 429 
Actually L,, leads to no violation of either C or P (1957). o 
invariance to any order. It is invariant under a charge- 6S. Glashow, Nuovo cimento (to be published). 
conjugation (or parity-reversing) operation which in- 'G. Feinberg and S. Weinberg, invited paper at 1960 
cludes a factor of i in the behavior of A relative to n. Winter Meeting of the American Physical Society at 
In this note, however, C and P are chosen with custom- Berkeley, California. 
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ANGULAR DISTRIBUTION OF LYMAN-a RA- 
DIATION EMITTED BY H (2S) ATOMS IN WEAK 
ELECTRIC FIELDS, William Lichten [Phys. 
Rev. Letters 6, 12 (1961). 


On page 12, line 7 from the bottom, read “In 
the present case J=3, m=+3,...” 
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In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


ELECTRICAL BREAKDOWN IN HYDROGEN AT 
LOW PRESSURES. A. L. Ward, Diamond 
Ordnance Fuze Laboratories, Washington, D. C., 
and Eifionydd Jones,* Department of Physics, 
University of Wales, University College of 
Swansea, United Kingdom (Received December 
13, 1960). 


Experimental and theoretical determinations 
of the static voltage-current characteristics, 
extending from the region of the Townsend (self- 
maintained) discharge to the normal glow dis- 
charge, have been carried out in hydrogen at 
low pressures (7-25 mm Hg). The calculations, 
made on an electronic computer, were based on 
the distortion of the electric field by space 
charge, and used the experimentally determined 
variation of both the primary and secondary 
Townsend ionization coefficients on the ratio 
of the field to pressure. Good agreement is ob- 
tained between the measured and calculated 
breakdown and glow voltages, and both the ex- 
perimental and theoretical curves of the charac- 
teristic are of similar shape. 


*Now at CERN, European Organization for Nuclear 
Research, Geneva, Switzerland. 


EXPERIMENTAL STUDY OF THE DIAMAG- 
NETISM OF GASEOUS PLASMAS WITH ELEC- 
TRON AND NUCLEAR SPIN RESONANCE TECH- 
NIQUES. T. C. Marshall and L. Goldstein, 
University of Illinois, Urbana, Illinois (Received 
December 15, 1960). 


It has been possible to observe the occurrence 
of diamagnetism in active discharges by meas- 
uring shifts in spin absorption resonance fre- 
quencies of foreign substances located near the 
plasma. Narrow spin electronic resonances 
in DPPH have been used to show that the mag- 
netization of a gaseous discharge in low-pressure 
mercury vapor increases linearly with applied 
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magnetic field up to about 40 gauss, where max- 
imization occurs. Observation of the phase 
change in the Larmor precession of protons in 

a strong homogeneous field showed that the 
diamagnetism in a modified Penning ionization 
gauge discharge increased linearly with power 
input and decreased approximately as 1/H. 
Irradiation of the plasma with high-power micro- 
waves with frequency near the free-electron 
gyrofrequency resulted in an increase of dia- 
magnetism of the discharge. The discharge 
magnetic moment ranged from -0.001 to -0.32 
erg/cm*-gauss, depending on experimental con- 
ditions. The sensitivity of the nuclear reso- 
nance technique is one part in 10°, and it per- 
mits observation of diamagnetism to be deferred 
until the discharge is over. The theory of plas- 
ma diamagnetism is summarized below. 


HIGHER RANDOM-PHASE APPROXIMATIONS 
IN THE MANY-BODY PROBLEM. H. Suhl and 
N. R. Werthamer, Bell Telephone Laboratories, 
Murray Hill, New Jersey, and Department of 
Physics, University of California, Berkeley, 
California (Received November 4, 1960). 


The usual random-phase approximation com- 
bined with an equations-of-motion technique for 
the many-electron problem is extended, yielding 
many of the known results of series summation 
methods in a straightforward manner. The 
method should apply to other types of many-body 
problems as well. 


REFLECTION OF SLOW ELECTRONS FROM 
TUNGSTEN SINGLE CRYSTALS, CLEAN AND 
WITH ADSORBED MONOLAYERS. P. Kisliuk, 
Bell Telephone Laboratories, Murray Hill, New 
Jersey (Received December 15, 1960). 


The reflection of electrons with kinetic energy 
up to a few electron volts from tungsten single- 
crystal surfaces is measured both on the clean 
surface and with adsorbed monolayers of nitro- 
gen and oxygen. For the clean surface, diffrac- 
tion from the lattice is responsible for a con- 
siderable part of the reflection in the thermionic 
range of energy. The magnitude of the reflection 
is such as to have a barely measurable effect on 
experimental tests of the thermionic emission 
equations. This technique permits continuous 
recording of the change in work function as gas 
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is adsorbed, yielding information about the 
kinetics of chemisorption and the surface dipoles 
due to the adsorbed gas atoms. 


ANHARMONIC ATTENUATION OF LOCALIZED 
LATTICE VIBRATIONS. P. G. Klemens, Wes- 
tinghouse Research Laboratories, Pittsburgh, 
Pennsylvania (Received December 9, 1960). 


Lattice modes localized about defects can in- 
terchange energy with the continuum of lattice 
waves by anharmonic interactions. The relaxa- 
tion time of a localized mode is calculated, 
taking account of cubic anharmonicities and 
using perturbation theory analogous to the treat- 
ment of three-phonon interactions. At zero 
temperature the relaxation time is typically of 
the order of 100 periods, but decreases with in- 
creasing temperature. 


POTENTIALS IN A CONDUCTOR OF VARYING 
CROSS SECTION. R. Jaggi, Research Labora- 

tory, International Business Machines Corpora- 
tion, Zurich, Switzerland (Received November 

21, 1960). 


A Bernoulli voltage, Vp, proportional to the 
square of the current J, may be expected in a 
conductor of varying cross section. Previous 
experiments to detect Vp are discussed. It is 
pointed out that in such experiments a Hall volt- 
age, Vzpxl 2, due to the magnetic field of the 
current J (“eigen- Hall effect” EHE), is super- 
posed on Vp. The ratio Vy/V zp is calculated. 
Experiments are performed on a bismuth sam- 
ple of varying cross section. The dependence 
of the measured voltages upon temperature and 
cross-sectional area shows that the EHE is 
dominant over the Bernoulli effect. 


EXCHANGE ANISOTROPY IN MIXED MAN- 
GANITES WITH THE HAUSMANNITE STRUC- 
TURE. I. S. Jacobs and J. S. Kouvel, General 
Electric Research Laboratory, Schenectady, 
New York (Received December 2, 1960). 


Magnetic properties are reported of powder 
compounds isomorphous to hausmannite, but 
containing partial to nearly complete substitution 
of diamagnetic Zn?* or Mg** ions in the tetra- 
hedral sites of Mn,O,. X-ray diffraction of 


these compounds reveals a single phase, whose 
structure corresponds to a tetragonally distorted 
spinel. Unidirectional anisotropy is present and 
is detected by the observation of hysteresis loops 
displaced along the field axis when the materials 
are cooled to low temperatures in magnetic 
fields of several kilo-oersteds. The unidirec- 
tional behavior is stable to reverse field pulses 
of 140 koe. An exchange anisotropy model is 
proposed involving interactions between ferri- 
magnetic and nearly antiferromagnetic regions 
brought about by the random distribution of the 
diamagnetic ions among the tetrahedral sites 
and the consequent magnetic inhomogeneity. 
This model is similar to that proposed for dis- 
ordered alloy systems. Magnetic viscosity ef- 
fects are observed which correspond in part to 
ordinary magnetic viscosity, but which are in- 
fluenced by the presence of unidirectional ani- 
sotropy. The identity of this anisotropy is not 
impaired by the ordinary viscosity effects. In 
addition, observations of certain differences in 
the hysteresis loops measured dynamically and 
quasi-statically reveal an “extraordinary” 
viscosity associated with a gradual breakdown 
of the unidirectional anisotropy. 


PARAMAGNETIC RESONANCE OF Gd°* IN 
ALO,. S. Geschwind and J. P. Remeika, Bell 
Telephone Laboratories, Murray Hill, New 
Jersey (Received December 7, 1960). 


The electron paramagnetic resonance spectrum 
of a small impurity of Gd** in Al,O, has been ex- 
amined at 24 kMc/sec. The over-all zero-field 
splitting of the ground state of 1.24 cm™ is the 
largest so far observed for Gd°*. The analysis 
of the spectra suggests that the Gd°*, which has 
twice the ionic radius of aluminum, essentially 
entered substitutionally for the Al°* (but dis- 
torted the environment in such a way as to ap- 
proach a condition of ninefold oxygen coordina - 
tion, the symmetry still remaining C,). The 
substitution of Gd°* for Al°* whose ionic radius 
is half as large would indicate that at impurity 
levels Jess than 0.02% matching of ionic radii is 
not an important criterion for incorporation into 
the lattice. 

Although there are two types of Al sites which 
are physically equivalent, the Gd°* entered 
selectively into one of these sites. This seem- 
ingly paradoxical result is ascribed to the dynam- 
ics of the crystal growth. The sites referred to 
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are actually inequivalent during the growth proc- 
ess as the Gd falls into place and are only equiv- 
alent in the grown crystal. 


CYCLOTRON RESONANCE IN GERMANIUM. 
Ralph R. Goodman,* Department of Physics, 
University of Michigan, Ann Arbor, Michigan 
(Received April 25, 1960; revised manuscript 
received January 9, 1961). 


Early cyclotron resonance for holes in Ge 
gave evidence of two resonant peaks associated 
with the so-called light and heavy holes. Luttin- 
ger and Kohn predicted, on the basis of a careful 
investigation of the theory of degenerate bands, 
that at low enough temperatures additional peaks 
should appear. Subsequent experiments con- 
firmed this prediction. In this paper a compari- 
son is made of the cyclotron resonance theory 
and the experiments of Fletcher, Yaeger, and 
Merritt. Values of the effective-mass constants 
which best fit the data are found. Resonant 
peaks additional to those found by Fletcher, 
Yaeger, and Merritt are predicted and discussed. 


*Now at Colorado State University, Fort Collins, 
Colorado. 


ELECTRON ENERGY BANDS OF ONE-DIMEN- 
SIONAL RANDOM ALLOYS. J. S. Faulkner* 
and J. Korringa, Department of Physics and 
Astronomy, The Ohio State University, Colum- 
bus, Ohio (Received December 5, 1960). 


A method for calculating the density of states 
for an infinite, one-dimensional random alloy 
is obtained by investigating the asymptotic be- 
havior of the trace of the “transmission” matrix 
which relates the values taken on by the wave 
function and its derivative at either end of the 
crystal. This matrix can be calculated if the 
potentials of the constituent A and B atoms, V A 
and Vg, are given. The equations are first 
derived for a very general case, and then the 
results of a calculation for an alloy in which the 
A and B atoms have equal concentrations is shown 
for the case that V4 and Vp are 6-function poten- 
tials. Certain generalizations of the method for 
treating other nonperiodic problems are dis- 
cussed briefly. 


*Now with the Department of Physics and the Quan- 
tum Theory Project, University of Florida, Gaines- 
ville, Florida. 
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NUCLEAR QUADRUPOLE RESONANCE IN AN 
ANTIFERROMAGNET. J. C. Burgiel,* V. Jac- 
carino, and A. L. Schawlow, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Re- 
ceived December 13, 1960). 


Nuclear quadrupole resonance techniques have 
been used to investigate magnetic interactions 
in the bromates and iodates of certain 3d transi- 
tion metal ions. In particular, for Ni(IO,).-2H,O 
there is an abrupt disappearance of the I’*"(+3/2 
+++1/2) transition at 3.08°K; at lower tempera- 
tures there is a large, temperature-dependent 
splitting. This behavior is attributed to the com- 
bined effects of an antiferromagnetic ordering of 
the Ni** electron spins and a hfs interaction of 
the nonlocalized Ni** spin magnetization with the 
I?” nuclear magnetic moment. Measurements of 
the pure quadrupole transition frequencies above, 
and of their temperature-dependent splittings 
below, the Néel temperature yield the quadrupole 
coupling constant, the asymmetry of the electric 
field gradient (EFG) tensor, the magnitude and 
orientation of the internal magnetic field relative 
to the principal axes of the EFG tensor, and the 
sublattice magnetization as a function of tem- 
perature. Qualitative experiments at 1.3°K indi- 
cate T,'*”«<100 sec and T,'*7~10~* sec. Meas- 
urements of the rf and dc magnetic susceptibility 
have been made and are consistent with these 
conclusions and, in addition, indicate the sudden 
onset of a spontaneous ferromagnetic moment as 
the temperature is lowered below T,- A possible 
isotope effect on the 7, of Ni(IO,),-2H,O was 
looked for but none was found. Cupric iodate and 
the bromates give no evidence of magnetic order- 
ing above 1.3°K. 


*Present address: Department of Physics, Massa- 
chusetts Institute of Technology, Cambridge, Massa- 
chusetts. 


| LATTICE THERMAL CONDUCTIVITY OF GER- 


MANIUM-SILICON ALLOY SINGLE CRYSTALS 
AT LOW TEMPERATURES. Arnold M. Toxen, 

International Business Machines Research Cen- 
ter, Yorktown, New York (Received December 

9, 1960). 


Thermal conductivity measurements are re- 
ported for five single-crystal Ge-Si specimens 
containing 0-7.56 atomic percent Si. The meas- 
urements were made under steady-state condi- 
tions and cover the temperature range 2-50°K. 
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The experimental results are compared to three 
theoretical models, those of Berman et al., 
Callaway, and Klemens; it is found that the data 
are best fit by Callaway’s model. Good agree- 
ment between experimental results and theoret- 
ical models is obtained by postulating only three 
sources of phonon scattering in the specimens: 
three-phonon processes, isotopic point-defect 
scattering by the germanium and silicon atoms, 
and boundary scattering. However, evidence is 
presented that boundary scattering occurs not 
only at the external surfaces of the specimens, 
but also at internal surfaces associated with 
microscale fluctuations of composition of the 
type reported by Goss, Benson, and Pfann. 


ELECTRICAL CONDUCTIVITY OF SINGLE- 
CRYSTAL CUPROUS OXIDE AT HIGH TEM- 
PERATURES. Robert S. Toth,* Rein Kilkson,T 
and Dan Trivich, Department of Physics and 
Department of Chemistry, Wayne State Univer- 
sity, Detroit, Michigan (Received November 14, 
1960). 


The electrical conductivity of single-crystal 
Cu,O0 was measured in the temperature range 
of 1100° to 500°C in oxygen pressures from 152 
mm to 10°° mm of Hg. The logo vs log P(O,) 
curves were found to be linear between the oxy- 
gen pressures of 50 mm and 10°? mm, with an 
average slope of 0.1420, or approximately 1/7. 
These curves exhibit a radical change in slope 
at O, pressures below 10°? mm. 

The plots of logo vs 1/T at constant oxygen 
pressure were found to be linear and the activa- 
tion energies obtained from the slopes of these 
plots have an average value of 0.65 ev at O, 
pressures between 50 mm and 10°? mm. At O, 
pressures of 10°* mm to 10°* mm, the activa- 
tion energy increases sharply to a value of 1.05 
ev. The activation energy obtained from the 
measurement of single-crystal Cu,O in air at 
temperatures from 1020°C to 1100°C was found 
to have an average value of 0.767 ev. 

An explanation for the physical significance of 
the activation energies obtained is suggested, 
and the models proposed to explain the depend- 
ence of the electrical conductivity on the O, 
pressure are considered. 


*Now at the Physics and Electronics Department, 
Ford Scientific Laboratory, Dearborn, Michigan. 

tNow at the Biophysics Department, Yale Univer- 
sity, New Haven, Connecticut. 








INFRARED CYCLOTRON RESONANCE IN InSb. 
E. D. Palik, G. S. Picus,* S. Teitler, and R. F. 
Wallis, U. S. Naval Research Laboratory, Wash- 
ington, D. C. (Received November 14, 1960; re- 
vised manuscript received January 23, 1961). 


Far-infrared cyclotron resonance absorption 
in n-type InSb has been measured to determine 
the variation of the conduction electron effective 
mass with magnetic field. At high magnetic 
fields the absorption was resolved into a strong 
line with a weaker satellite line at lower photon 
energy and a broad, weak absorption at still 
lower energy. The interpretation of this struc- 
ture in terms of the conduction band properties 
of InSb is discussed. 


*Now at Hughes Semiconductors Products Division, 
Newport Beach, California. 


NECESSITY AND EXPERIMENTAL CONSIS- 
TENCY OF ANTIFERROMAGNETIC GROUND 
STATE WITHOUT LONG-RANGE ORDER. G. W. 
Pratt, Jr., Lincoln Laboratory, Massachusetts 
Institute of Technology, Lexington, Massachu- 
setts (Received December 5, 1960). 


A proof is given that long-range order meas- 
ured as the average value of the z component of 
sublattice magnetization of an antiferromagnet 
must be zero in the exact ground state if that 
state is nondegenerate and the Hamiltonian is 
invariant under time reversal. The neutron dif- 
fraction magnetic cross section is shown not 
to depend on such a measure of long-range order. 


PRESSURE EFFECT ON SUPERCONDUCTING 
LEAD. M. Garfinkel* and D. E. Mapother, De- 
partment of Physics, University of Illinois, 
Urbana, Illinois (Received November 29, 1960). 


Techniques are described for measuring the 
effect of hydrostatic pressure on the critical 
field, H,, of superconducting Pb. Pressures up 
to 650 atm were applied using solid helium as 
the pressure fluid. Observations were made 
from about 7 to 1°K, and values of dH,/dP, 
dT,./dP, and the temperature variation of 
(8H ./eP)7 are reported. From these data the 
value of (1/y*)(dy*/dP) is deduced, where y* 
is the temperature coefficient per unit volume 
of the normal electronic specific heat. The ob- 
served data are accurately represented over the 
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full range of measurement by the equation 

HP, T)=H,(PY (t), where t=T/T, and f(t) is 
independent of pressure. The “similarity prin- 
ciple” requirement, H,(P)/T,(P)=const, is 
shown to be invalid for Pb. The results provide 
the basis for a discussion of the pressure effects 
on the net interaction potential, V, of the Bar- 
deen, Cooper, Schrieffer theory and the density 
of electronic states near the Fermi surface. 


*Present address: General Electric Research 
Laboratory, Schenectady, New York. 


GALVANOMAGNETIC PROPERTIES OF n-TYPE 
CdAs,. A. S. Fischler, International Business 
Machines Corporation, Poughkeepsie, New York 
(Received November 21, 1960). 


Galvanomagnetic measurements on oriented, 
single crystals of n-type CdAs,, a noncubic 
semiconductor, indicate the surfaces of constant 
energy to be ellipsoids of revolution, located 
along the symmetry axis of the crystal system. 
The ratio of electronic mobility is found to be 
by/by ~4, from Hall and resistivity data. Mag- 
netoresistance measurements confirm this con- 
duction band model and indicate the scattering 
to be due primarily to acoustical lattice modes 
with some degree of impurity scattering. 


PHOTOTHERMAL EFFECT IN SEMICONDUC- 
TORS. W. W. Gartner, C. B. S. Laboratories, 
Stamford, Connecticut (Received November 17, 
1960; revised manuscript received December 27, 
1960). 


When a sample of semiconducting material is 
illuminated, pairs of excess carriers are gener- 
ated which diffuse through the material according 
to the density gradients established. Each pair 
carries an energy approximately equal to the 
band gap of the material. This energy is depos- 
ited where the excess electron recombines with 
a hole and causes local heating of the lattice. A 
temperature distribution will therefore be estab- 
lished in the sample which depends on the charac- 
teristics of optical absorption and bulk and sur- 
face recombination in and on the sample. This 
establishment of a temperature distribution in a 
solid by optically excited diffusing and recombin- 
ing carriers is called the photothermal effect. 
The paper gives a formulation of the theory gov- 
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erning the photothermal effect, and the case of 
small temperature elevations in an infinite slab 
is worked out in detail. 


STATISTICAL MECHANICS OF FERROMAG- 
NETISM; SPHERICAL MODEL AS HIGH- DEN- 
SITY LIMIT. R. Brout, Laboratory of Atomic 
and Solid-State Physics and Department of Phys- 
ics, Cornell University, Ithaca, New York (Re- 
ceived November 15, 1960). 


It is shown that there is a class of graphs of 
the Ising model (or Heisenberg model for 
T> Tcurie) which is comprised of cycle graphs 
plus some excluded-volume effects which sum to 
the spherical model. The spherical model, suit- 
ably generalized for T< Tcurie’ was conjectured 
in a previous work to be the high-density limit 
of the Ising model, correct to 1/z. z“° meas- 
ures the range of the exchange potential. This 
is now proved by examining the omitted graphs. 
The error is shqwn to be O(1/z?). 


MAGNETIC TORQUE CURVES FOR A SINGLE 
CRYSTAL OF THULIUM ORTHOFERRITE 
(TmFeO,). Chinji Juroda, Tomonao Miyadai, 
Akira Naemura, Nobukazu Niizeki, and Hisao 
Takata, Electrical Communication Laboratory, 
Nippon Telegraph and Telephone Public Corpora- 
tion, Musashino-shi, Tokyo, Japan (Received 
November 21, 1960). 


The magnetic anisotropy of a thulium ortho- 
ferrite single crystal was investigated by means 
of an automatically recording torque-magneto- 
meter in the range between 77 and 300°K. It was 
found that at room temperature the direction of 
weak magnetization, 0,, associated with an anti- 
ferromagnetic canted spin arrangement is par- 
allel to the c axis of orthorhombic structure, 
and that at about 100°K a magnetic transition 
occurs, where the direction of 0, changes from 
the c axis to the a axis. The transition tempera- 
ture obtained is the highest among those re- 
ported before on other orthoferrites. The tor- 
que curves obtained above the transition tem- 
perature are similar to those by Sherwood et al. 
and can be explained by their model. The value 
of o, is equal to 0.06 » p/mole at room tempera- 
ture. At 77°K below the transition temperature, 
however, inexplicable torque curves were ob- 
tained in the b plane. The curves showed a 
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small jump at the c axis superimposed on a 
sin26 curve with large amplitude. The jump 
amplitude is almost independent of the applied 
field in the range of 8500 to 3000 oe. This be- 
havior cannot be explained by Sherwood’s model, 
but it will give information on the interaction 
between thulium and ferric ions. 


GALVANOMAGNETIC EFFECTS IN n-Ge IN THE 
IMPURITY CONDUCTION RANGE. R. J. Sladek 
and R. W. Keyes,* Westinghouse Research La- 
boratories, Pittsburgh, Pennsylvania (Received 
December 13, 1960). 


Measurements of the magnetoresistance and 
magnetic field dependence of the Hall coefficient 
of several samples of n-type germanium in the 
impurity conduction range have been made em- 
ploying magnetic field strengths up to 28 kgauss. 
The magnitude and the crystalline anisotropy of 
the magnetoresistance are interpreted in terms 
of the changes in the donor wave functions which 
are produced by the magnetic field. The field 
dependence of the Hall coefficient is interpreted 
as a magnetoresistance effect of the conduction 
band. 


*Present address: IBM Research Laboratory, 
Poughkeepsie, New York. 


EFFECT OF CONFIGURATION MIXING AND 
COVALENCY ON THE ENERGY SPECTRUM OF 
RUBY. S. Sugano* and M. Peter, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Re- 
ceived December 1, 1960). 


For the purpose of improving the analysis 
hitherto done of the optical and microwave spec- 
trum of ruby, a calculation has been performed 
in the strong cubic field scheme, taking into 
account the effect of configuration mixing of the 
higher excited t,’e states into the ¢,° states. In 
the calculation, covalency of the ¢, and e elec- 
trons is also introduced in a simplified fashion 
besides the spin-orbit interaction, trigonal field, 
and Zeeman energy. 

The result shows that configuration mixing and 
covalency play a very important role in giving 
zero-field splittings and g values of the ¢,° states. 
It is also found that there is not much difference 
between the degrees of covalency for the ¢, and 
e electrons, although they are fairly large for 
both electrons. The best zero-field splitting of 











the ground state thus obtained is 0.24 cm™ with 
the correct sign. 


*On leave from the Department of Physics, Tokyo 
University, Tokyo, Japan. 


SUPERCONDUCTING TRANSITION TEMPERA- 
TURE OF HIGH-PURITY TANTALUM METAL. 
J. G. C. Milne, Department of Physics, The 
University of Sheffield, Sheffield, England (Re- 
ceived October 3, 1960). 


The superconducting transition temperature 
of a high-purity sample of tantalum has been 
measured and a zero-field value of 4.457°K 
+ 0.003°K obtained. The ratio of room tempera- 
ture resistivity p,,, to the residual resistivity 
P4.2 gave a measure of the purity of the sample, 
a value of 250 being obtained. These results 
compare very well with those of Budnick. The 
value of (dH,./dT) near the zero-field transition 
temperature was -358 oersteds/°K. An approxi- 
mate theoretical estimate similar to that done 
by Pippard for tin shows that small amounts of 
impurity as measured by the residual resistance 
ratio will give a decrease in transition tempera- 
ture of about half that observed experimentally. 


PHOTOCONDUCTIVITY AND TRAPPING IN 
SILVER CHLORIDE CRYSTALS. Albert M. 
Gordon,* Laboratory of Atomic and Solid State 
Physics and Department of Engineering Physics, 
Cornell University, Ithaca, New York (Received 
December 15, 1960). 


Photoconductivity and trapping were investi- 
gated in pure, cuprous- and nickel-chloride- 
doped, and darkened silver chloride crystals. 
The photoconductivity was measured primarily 
at 88°K. The initial photoresponse increases 
with rising absorption constant, peaks at wave- 
lengths for which the absorption constant was 3 
to 5cm-~', and falls rapidly to 1/20 to 1/10 of 
the peak value at wavelengths at which the ab- 
sorption constant was 50 to 100 cm™'. This de- 
crease in photoresponse at short wavelengths is 
explained in terms of trapping and recombination 
through centers in a surface region. The effect 
of irradiation and the filling of traps was investi- 
gated. Values of the schubweg for electrons in 
the bulk material were obtained from these ex- 
periments. The cuprous-chloride-doped samples 
had a long-wavelength tail on the photoresponse 
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curve corresponding to the long-wavelength tail 
observed on the optical absorption. 

The photoconductivity does not show any peaks 
in the infrared. There was, however, a photo- 
response with a threshold at approximately 
10000 A in samples irradiated at 88°K anda 
peak for those irradiated at room temperature. 
This response was attributed to photoemission 
from free silver introduced unavoidably during 
the preparation of the samples. The warming 
of samples irradiated at 88°K to fill traps pro- 
duced no measurable thermoluminescence. 
Electrical warming curves (measurement of cur- 
rent as a function of time during warming) dis- 
closed a number of trapping levels. All samples 
except the nickel-doped sample showed peaks at 
115°K, 140°K, and 180°K related to activation 
energies of approximately 0.20 ev, 0.28 ev, and 
0.45 ev. The nickel-doped sample showed only 
peaks at 115°K and 180°K. The cuprous- doped 
samples and the darkened pure sample showed 
a peak at 160°K corresponding to a thermal activa- 
tion energy of approximately 0.36 ev. On some 
samples another peak was visible at 240°K with 
a thermal activation energy of approximately 
0.62 ev. Possible interpretations for these vari- 
ous peaks are discussed. 


*Present address: Department of Physiology, Uni- 
versity College, London, England. 


TEMPERATURE DEPENDENCE OF THE WORK 
FUNCTION OF METALS (Mo, Ni). George Comga 
Adrian Gelberg, and Beatrice Iosifescu, Institute 
of Atomic Physics, Bucharest, Rumania (Re- 
ceived August 29, 1960; revised manuscript re- 
ceived December 15, 1960). 


The temperature coefficient of the work func- 
tion (TCWF) of polycrystalline Mo and Ni has 
been measured in the temperature ranges 600- 
1100°K and 475-1025°K, respectively, using the 
electron beam method. Special precautions have 
been taken in order to avoid systematic errors 
due to the effect of residual gases, of stray mag- 
netic fields, of incorrect temperature measure- 
ment, of cathode temperature variations, etc. 
During determinations made upon Ni surfaces 
the terrestrial magnetic field has been compen- 
sated. Random error sources have been likewise 
reduced as much as possible. The main results 
obtained on surfaces outgassed for many thou- 
sands of hours at pressures of 1x10-*°- 3x107® 
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torr are 
(yy,./47) ay 7 (7-84+ 0.07) x 10°* ev/°K, 


={- -5 ° 
(@d,/4T) 7, Seis’ (-3.12 + 0.05) x 10-5 ev/°K. 


At the Curie point of Ni, a theoretically ex- 
pected variation of the TCWF has been observed: 


I/O pcg ~ Y/4 750 av 


=(-0.99+ 0.17) x 10-5 ev/°K, 


but no jump of the work function has been found 
at this point. Special experiments have shown 
that the results have not been significantly af- 
fected by residual gases. A relation between the 
TCWF of the various facets of these surfaces has 
been obtained. 


SPECIFIC HEAT OF YTTRIUM IRON GARNET 
FROM 1.5°K TO 4.2°K. Samuel S. Shinozaki,* 
Department of Physics, University of California, 
Berkeley, California (Received December 19, 
1960). 


The specific heats of two samples of yttrium 
iron garnet have been measured between 1.5 and 
4.2°K. The data have been analyzed into lattice- 
and spin-wave contributions characterized, re- 
spectively, by the Debye temperatures 0, = 538°K, 
©,=567°K, and by D, =0.81x10-** erg cm’, D, 
= 0.85x10-** erg cm’, where D is defined by the 
dispersion relation for spin waves, hw = Dk?. 


*Present address: Ford Scientific Laboratory, 
Dearborn, Michigan. 


ELECTRON CAPTURE FROM He(1s?) BY PRO- 
TONS. Robert A. Mapleton, Geophysics Re- 
search Directorate, Air Force Cambridge Re- 
search Laboratories, Air Research and Develop- 
ment Command, L. G. Hanscom Field, Bedford, 
Massachusetts (Received October 20, 1960). 


The two equivalent forms of Born’s approxi- 
mation, prior and post, are used to calculate 
the electron capture cross section for protons 
incident on He(1s?). These cross sections are 
calculated for capture into eleven different final 
states in the energy range 12.5 kev to 1 Mev. 
Although a rather crude wave function, (Z*/na,°) 
x exp[-(Z/a,)(7,+7,)] (Z = 1.6875), is used for 
He, the prior and post total capture cross sec- 
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tions do not differ by more than 20% over the 
energy range investigated. Estimates of the 
sum of the cross sections for capture into all 

s states of the hydrogen atom for the two resid- 
ual ions, He*(1s) and He*(2s), are obtained 
from an adaption of the s-state sum rule as 
given in the paper of Jackson and Schiff. As in 
this work of Jackson and Schiff, it is found that 
the s states provide the major contribution to 
the total capture cross section. The calculated 
cross sections agree fairly well with the ex- 
perimental values. The cross section for cap- 
ture into the state Het(1s)+H(1s) is roughly 
2.5 times larger than the values obtained by 
Bransden, Dalgarno, and King. 


IONIZATION PRODUCED BY ATOMIC COLLI- 
SIONS AT kev ENERGIES. Il. J. B. Bulman* 
and A. Russek, University of Connecticut, Storrs, 
Connecticut (Received December 13, 1960). 


The electron evaporation model of the collision- 
ionization process that occurs when atoms collide 
at high energies is extended to include atoms 
with from two to eight electrons in the outer 
shell. Application of the model to data from col- 


lisions of N* on Ar and Ne* on Ar gives evidence 
for a resonant electron capture effect taking 
place in high-energy violent collisions which 

was heretofore masked by the multiple ioniza- 
tion consequent on such collisions. 


*Now at Central Connecticut State College, New 
Britain, Connecticut. 


STRUCTURE IN THE IONIZATION NEAR THRESH- 
OLD OF RARE GASES BY ELECTRON IMPACT. 

S. N. Foner and B. H. Nall, Applied Physics La- 
boratory, The Johns Hopkins University, Silver 
Spring, Maryland (Received August 1, 1960; re- 
vised manuscript received January 25, 1961). 


Ionization efficiency curves for xenon, krypton, 
and argon have been studied with an electron en- 
ergy analyzer. The electron energy distribution 
was measured and the absolute voltage scale de- 
termined in each experiment. The results of 
these studies (1) favor a linear threshold ioniza- 
tion law over a 1.127 power law, and (2) show 
that the data cannot be explained simply by ioni- 
zation processes with onsets at the *P,, and 
*P. ground states of the ion, but can be well 
fitted by a series of linear processes. The ioni- 


zation potentials obtained by extrapolating ac- 
cording to a linear threshold law agree with 
spectroscopic values to within 0.02 ev. New on- 
sets in argon were observed at about 0.64 v and 
1.27 v above threshold. The observed structures 
in the rare gases are not readily explained by 
auto-ionization and no alternative explanation is 
offered. The structures observed in these experi- 
ments are compared with the results obtained by 
other “high-resolution” techniques. This com- 
parison is complicated by the disparity in the 
published data on onset energies, and by the 
even greater disagreement on the relative pro- 
babilities for the various ionization processes. 
An independent check on consistency of data was 
made by comparison with “low-resolution” data 
obtained on a conventional mass spectrometer. 
The present data are in excellent agreement with 
the lower resolution data, while some of the 
other “high-resolution” data are not. 


PROTON-HYDROGEN SCATTERING SYSTEM. 
Marvin H. Mittleman, Lawrence Radiation Labor - 
atory, University of California, Livermore, 
California (Received August 26, 1960). 


The impact parameter treatment of the scatter- 
ing of protons by hydrogen is derived and is 
shown to be valid for energies greater than a few 
volts. A novel treatment of the resultant equa- 
tions is given which significantly modifies pre- 
viously obtained results for inelastic scattering 
and charge-exchange scattering. 


SPECTRA INDUCED BY 200-kev PROTON IM- 
PACT ON HELIUM. R. H. Hughes, R. C. War- 
ing,* and C. Y. Fan,t Physics Department, Uni- 
versity of Arkansas, Fayetteville, Arkansas 
(Received November 3, 1960; revised manu- 
script received January 9, 1961). 


Spectra induced by 200-kev proton impact on 
helium have been observed in the spectral re- 
gion of 43500 A to ,6000 A. 'S states appear to 
be strongly excited. Absolute cross sections for 
the direct excitation of the 4 'S and 5 'S states 
of neutral helium were determined as well as 
the simultaneous ionization and excitation cross 
section for helium into the n=4 state of Het. Of 
the more intense lines only the 2 'P-n 'S lines 
and the Hell 4686 line behaved linearly with 
pressure within experimental error. Triplet 
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spectra were observed in which the dominant 
feature was the 2 *P-n °D lines. The popula- 
tions of the 4 °D and 4 'D states, in particular, 
were analyzed as a function of a direct mecha- 
nism and collision of the second kind which seem 
to fit the data fairly well. A very weak Doppler- 
shifted H, line was detected. If this is inter- 
preted to be produced by charge exchange, then 
the cross section for electron capture into the 
n=4 state of hydrogen is estimated to be of the 
order of 8x10°™** cm’. 


*Now at the University of Kansas City, Kansas City, 
Missouri. 

tNow at the Laboratories for Applied Science, Uni- 
versity of Chicago, Chicago, Illinois. 


VARIATION- PERTURBATION METHOD FOR 
EXCITED STATES. Oktay Sinanoglu,* Depart- 
ment of Chemistry and Lawrence Radiation Labo- 
ratory, University of California, Berkeley, 
California (Received July 22, 1960). 


The first-order wave function, X,, in the per- 
turbation method can be obtained by a variational 
principle instead of summing the usual infinite 
series with a large continuum contribution. For 
a ground state or the lowest state of a given 
symmetry suitable trial functions, x, are chosen 
to attain E,, the second-order contribution to 
the energy, as a minimum. This method is ex- 
tended here to any excited state, m. To obtain 
x,”™, the expression 


m.~ m 


~ m m 
E, -{2(2, ,H,-E, )X, ) 


~ m m.~ m m 
+(X, »H,-E, )X, Y 2B, 


is to be minimized with X,™ in the form 
k,_k m 
m-1%) (%) ,H,®% ?) 


x =x," + p>» 
k=1 





m > 
“~— 
and X,”” orthogonal to the known unperturbed 
functions of the states lower than m. The X, 
gives also the third-order energy. 


*Present address: Sterling Chemistry Laboratory, 
Yale University, New Haven, Connecticut. 
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PERTURBATION THEORY OF MANY-ELECTRON 
ATOMS AND MOLECULES. Oktay Sinanoglu,* 
Department of Chemistry and Lawrence Radia- 
tion Laboratory, University of California, 
Berkeley, California (Received July 25, 1960). 


Perturbation theory with operator techniques is 
applied to a nondegenerate many-electron system 
taking the electron repulsions, ));% j rg}, as 
the perturbation. The first-order wave function, 
X,, is obtained rigorously in terms of the first- 
order wave functions of independent two-electron 
systems. The pair functions containing nuclear 
parameters can be solved for individually by 
variational or other methods, then used in vari- 
ous atoms or molecules. For example, the Li 
atom is built up completely from the (1s)? 'S, 
(1s2s) *S and °S states of Lit. The X, gives the 
energy to third order and as an upper limit to 
the exact E. The E, is equal to the sum of com- 
plete pair interactions plus many-body terms of 
two types: (a) “cross polarization” which exists 
even in no-exchange intermolecular forces, and 
(b) Fermi correlations. 


*Present address: Sterling Chemistry Laboratory, 
Yale University, New Haven, Connecticut. 


APPLICATION OF THE IMPULSE APPROXIMA- 
TION TO THE SCATTERING OF ELECTRONS 
BY ATOMS: I. INELASTIC SCATTERING BY 
HYDROGEN ATOMS. R. Akerib, Convair, San 
Diego, California, and S. Borowitz, Physics 
Department, New York University, New York, 
New York (Received November 18, 1960). 


A calculation of the ionization cross section 
and the excitation cross section to the 2S and 
2P states of hydrogen atoms by collision with 
electrons is carried through by the use of an 
impulse approximation. The results are then 
compared to the experimental data and to var- 
ious other theoretical calculations. The cross 
sections obtained compare favorably with ex- 
periment. The calculations carried out by these 
methods are no more complicated than the usual 
approximation methods and are easily adaptable 
for use with more complicated atoms. 


ELASTIC SCATTERING AND REACTIONS OF 
PROTONS ON O”*. R. R. Carison, C. C. Kim,* 
J. A. Jacobs,t and A. C. L. Barnard,} Depart- 
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ment of Physics and Astronomy, State Univer- 
sity of Iowa, Iowa City, Iowa (Received Decem- 
ber 14, 1960). 


The elastic scattering 0'*(p,p)O** and the re- 
actions O**(p,p’y)O**, O**(p, a,)N**, and 
0(P, a; 2¥1,2)N** were studied using a thin gas 
target with the State University of Iowa electro- 
static generator. Absolute differential cross 
sections were measured for the two laboratory 
angles 86.8° and 159.5° in the incident proton 
energy range 790 to 3550 kev, and angular distri- 
butions for a, and p were measured at several 
energies. Relative yield curves of gamma rays 
were obtained in the same energy range as 
above. Two F'® levels were observed which 
have not been previously reported and some new 
decay modes for previously known levels were 
observed. From consideration of the detailed 
shape of the elastic-scattering anomalies and 
the angular distributions, spin and parity as- 
signments were made to some F’® levels. 


*Now at Physics Department, University of Southern 
California, Los Angeles, California. 

TNow at Physics Department, Virginia Polytechnic 
Institute, Blacksburg, Virginia. 

tNow at Physics Department, Rice University, 
Houston, Texas. 


GAMMA-GAMMA DIRECTIONAL CORRELA- 
TIONS IN Cs***, Atam P. Arya,* Department of 
Physics, The Pennsylvania State University, 
University Park, Pennsylvania (Received Novem- 
ber 23, 1960). 


Directional correlation measurements have 
been made on the 356—82 kev, 301-82 kev, and 
80-82 kev gamma-gamma cascades in Cs*** 
following the decay of 8-year half-life Ba*** with 
a coincidence scintillation spectrometer using 
Nal detectors. The observed correlation func- 
tions are: 


W(0@) =1+(0.042 + 0.005)P,(cosé) 
~(0.0041 + 0.0038)P,(cos@), 
W(@) =1 -(0.0257 + 0.011)P,(cosé@) 
-(0.0002 + 0.008)P,(cos@), 
W(6) =1+(0.0487 + 0.017)P,(cosé) 
+(0.0011 + 0.012)P,(cosé@), 


respectively, for the three cascades. These 
gamma-gamma directional correlations were 





found to be consistent with spin assignments of 
7/2*, 5/2*, 3/2+, 3/2*, and 1/2* to the levels 
at the ground state, 82 kev, 162 kev, 383 kev, 
and 438 kev in Cs***. The probable assignment 
of the multipolarities for different gamma rays 
is as follows: The 356-kev gamma ray is pure 
E2. The 82-kev gamma ray is a mixture of 
(96.54 0.5)% M1 and (3.54 0.5)% £2 with 6,, 
=-0.190+0.014. The 301-kev gamma ray can 
have one of two possible mixtures: either 6,,, 
=+0.123 + 0.004 with a mixture of (98.5+1.0)% 
M1 and (1.5+1.0)% £2, or 6,,, =3.977+1.02 with 
a mixture of (6+2)% M1 and (9442)% £2. The 
value of 6,,, = -3.977+ 1.02 is more probable. 
The 80-kev gamma ray is also found to have 
two possible values of 6,,: either 6,,=+0.468 

+ 0.094 with a mixture of (824 6)% M1 and 

(18+ 6)% £2, or 6,,=+7.0 with a mixture of 
(2.0+1.5)% M1 and (98.04 1.5)% £2. 


*Permanent home address: 3623, Chawri Bazar, 
Delhi-6, India. 


SIMPLE REALISTIC TREATMENT OF NUCLEAR 
DIRECT-INTERACTION PROCESSES. Ian B. 
McCarthy* and Derek L. Pursey,t Department 

of Physics, University of California, Los 
Angeles, California (Received October 27, 1960; 
revised manuscript received January 24, 1961). 


Physical arguments are used to predict quali- 
tatively the effect on direct-interaction differ - 
ential cross sections of the distortion of the 
wave functions of the scattered particle. These 
qualitative predictions are confirmed by calcula- 
tions using a simple but fairly realistic model 
for the wave-function distortion in (a, a’) scat- 
tering. The model used is based on examination 
of the properties of optical-model wave functions. 
Good fits to experimental data are found using 
the model for (a, a) scattering in the energy 
range 20-40 Mev for scattering angles less than 
90°. Features of direct-interaction processes 
involving nucleons are interpreted in terms of a 
focus in the optical-model wave functions for 
these particles, but detailed calculations are 
not presented. 


*Permanent address: Department of Mathematical 
Physics, University of Adelaide, Adelaide, S. Austra- 
lia. 

Ton leave of absence from University of Edinburgh, 
Edinburgh, Scotland. Present address: Institute for 
Advanced Study, Princeton, New Jersey. 
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DECAY OF K® and Sc“. John D. McCullen* and 
J. J. Kraushaar, University of Colorado, Boulder, 
Colorado (Received December 2, 1960). 


The low excited states of Ca** and Ca* have 
been studied in the decay of K** and Sc**, with 
special emphasis on the observation of weakly 
populated states. In the K** decay, gamma rays 
have been seen at 0.31 Mev (1.1%), 0.49 Mev 
(<0.1%), 0.60 Mev (0.1%), 0.90 Mev (0.1%), 
1.02 Mev (0.1%), 1.52 Mev (100%), 1.92 Mev 
(0.3%), and 2.42 Mev (0.2%). The coincidence 
sequence of the transitions has been measured 
and a level scheme constructed. 

In the Sc** decay, gamma rays have been seen 
with energies and intensities of 0.68 Mev (3.2%), 
1.02 Mev (3.1%), 1.12 Mev (4.7%), 1.16 Mev 
(100%), 1.50 Mev (1.7%), 1.72 Mev (0.8%), 2.28 
Mev (0.2%), and 2.69 Mev (0.2%). Coincidence 
measurements were also taken for this isotope 
to clarify the cascade sequences, and a level 
scheme was constructed. A search was made 
for low-energy conversion electrons (E < 1 Mev) 
in an effort to establish the existence or non- 
existence of a low lying 0* state in Ca**, whose 
analog occurs as the second excited 1.84- Mev 
state in Ca**. No such conversion electrons 
were seen, either by electron spectrometer 
studies or by electron-delayed gamma-ray coin- 
cidence measurements. An upper limit of 0.05% 
of the total decay of Sc“ was put on the popula- 
tion of such a state. 


*Now at Princeton University, Princeton, New 
Jersey. 


INELASTIC SCATTERING OF DEUTERONS 
FROM THE MAGNESIUM ISOTOPES. A. G. 
Blair* and E. W. Hamburger, University of 
Pittsburgh, Pittsburgh, Pennsylvania (Received 
December 2, 1960). 


Natural and enriched magnesium targets were 
bombarded with 15-Mev deuterons from the 
University of Pittsburgh cyclotron. The reaction 
products were magnetically analyzed and detected 
by a scintillator or by nuclear emulsions. Angu- 
lar distributions were obtained for the deuterons 
inelastically scattered from the Mg” 1.37-, Mg” 
1.61-, and Mg”® 1.83-Mev states. The results 
are compared to the predictions of the plane- 
wave Born-approximation theory and of the in- 
elastic diffraction-scattering model. The curves 
obtained from the Born approximation give better 
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over-all correspondence with the experimental 
points. The inelastic diffraction-scattering 
model, however, allows one to extract directly 
the effective values of the nuclear deformation 
parameter 8. One obtains |8| =0.20, 0.19, and 
0.17 for Mg**, Mg**, and Mg”, respectively. 
The spectra of deuterons inelastically scattered 
from Mg” and Mg” were also observed, at 6), 
= 12°, 30°, and 60°. The only large cross sections 
in Mg”® were those for the 1.61-Mev (7/2*) level 
and for a level near 3.4 Mev. The strength of 
the reaction to the latter level suggests that it is 
the 9/2* member of the ground-state rotational 
band which, in analogy with Al*®, should appear 
at approximately this energy. The results tend 
to confirm the selection rule that favors collec- 
tive excitations over single-particle excitation 
in inelastic scattering. In Mg” strong scattering 
was observed only from the first two excited 
states. A previously unreported Mg” state was 
found at the excitation energy of 3.614+ 0.020 Mev. 


*Now at Los Alamos Scientific Laboratory, Los 
Alamos, New Mexico. 

Now at Faculdade de Filosofia, Ciéncias e Letras 
(Departamento de Fisica), Universidade de Sao Paulo, 
Si Paulo, Brazil. 


NUCLEAR SPIN OF Ho’. wW. J. Childs and L. S. 
Goodman,* Argonne National Laboratory, Argonne, 
Illinois (Received December 2, 1960). 


The hyperfine structure of Ho’ has been ex- 
amined by means of the atomic-beam magnetic- 
resonance technique. The atomic ground state 
appears to be ‘/,,., and g (y has the measured 
value 1.19509+ 0.00007, in close agreement with 
the Russell-Saunders value of 6/5. Only one re- 
sonance is observed, and its transition frequency 
for 0.15 <H < 150 gauss is proportional to the 
magnetic field strength to within experimental 
error. The simplest interpretation is that the 
nuclear spin is J=0. It is shown that if J=1, the 
magnetic hyperfine interaction constant a must 
be less than 5 kc/sec. 


*Currently on leave as a Guggenheim Fellow at the 
University of Heidelberg, Heidelberg, Germany. 


(p,d) PICKUP REACTIONS IN LIGHT NUCLEI. 
E. F. Bennett,* Palmer Physical Laboratory, 
Princeton University, Princeton, New Jersey 
(Received June 20, 1960). 


Using a thin proportional counter as a velocity 
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selector in conjunction with a Nal crystal to 
measure energy, deuterons from proton- induced 
reactions in some light nuclei were studied. The 
detection system was capable of presenting an 
essentially undistorted spectrum of deuterons 

in the presence of protons of the same energy 
and considerably more intense. Angular distribu 
tions of deuterons from C'5(p,d)C'*, N"*(p,d)N’, 
N‘*(p,d)N*, F'°(p,d)F"*®, Mg*(p,d)Mg™*, and 
P*!(p,d)P*° were taken. Butler curves were cal- 
culated to fit the experimental distributions and 
level widths extracted. 


*Now at Argonne National Laboratory, Argonne, 
Illinois. 


CARBON-13 NEUTRON TOTAL CROSS SECTION. 


H. O. Cohn, J. K. Bair, and H. B. Willard, Oak 
Ridge National Laboratory, Oak Ridge, Tennes- 
see (Received December 9, 1960). 


The total neutron cross section of carbon-13 
was measured for neutron energies from 110 
kev to 9 Mev and from 16 Mev to 23 Mev. Four 
narrow resonances were observed as well as 
broad resonance structure above 3 Mev. The 
most probable spin assignments based on the 
resonance heights and widths are as follows: 


= 153+ 5kev, [=13 kev, J=1*; 


Ep =17514 8kev, I =20 kev, 


Ep =2432+10 kev, [=17 kev, 


Ep =2454+10 kev, F=10 kev, 


NUCLEAR INTERACTIONS AND MEAN FREE 
PATHS OF PROTONS, NEUTRONS, AND a 
PARTICLES AT ENERGIES AROUND 250 Bev/ 
NUCLEON. E. Lohrmann, M. W. Teucher,* and 
Marcel Schein,! Department of Physics, Univer- 
sity of Chicago, Chicago, Illinois (Received 
November 2, 1960). 


Nuclear interactions of protons, neutrons, a 
particles, and heavier nuclei of average energy 
250 Bev/nucleon were studied in nuclear emul- 
sion. The sources of these particles were frag- 
mentations of heavy primary nuclei of the cosmic 
radiation. Their energy was determined from 
multiple scattering measurements. The interac- 


tion mean free path for protons is 41+ 10 cm, 

for a particles 27+ 7 cm. The mean free path 
shows no significant change compared with meas- 
urements at lower energies. The mean number 
of shower particles, (Ns); depends appreciably 
on the mass of the target nucleus. Our best 
estimate for nucleon-nucleon collisions at 250 
Bev is (N,)=8.841.9. A detailed comparison of 
the estimate of the primary energy obtained 
from the angular distribution of shower particles 
with the true primary energy is carried out. 

The angular distribution of the shower particles 
will, in an individual case, give a quite unre- 
liable value for the primary energy. On the 
average, the angular distribution method will 
overestimate the true primary energy by a fac- 
tor of 1.3 for interactions with a number of heavy 
prongs N;, <5. If N,>5, the angular distribution 
will underestimate the true energy in the average 
by a factor of 1.8. The angular distributions can 
be transformed into a system in which they are 
symmetric. This is even true for collisions with 
heavy target nuclei (N;,>5). The results for 
alpha-particle and heavy nucleus collisions are 
quite similar. 

The inelasticity for the proton and neutron in- 
teractions shows large fluctuations for individual 
events. It depends weakly on the number of 
shower particles and on the mass of the target 
nucleus. Its mean value is 50%. The mean 
value for the alpha-particle collisions is 22 %. 


*Now at the Physikalisches Staatsinstitut, Univer- 
sity of Hamburg, Hamburg, Germany. 
TDeceased. 


NUCLEAR INTERACTIONS OF PROTONS, 
NEUTRONS, AND SHOWER PARTICLES OF 
VERY HIGH ENERGY IN NUCLEAR EMULSION. 
A. Barkow,* B. Chamany,! D. M. Haskin, P. L. 
Jain,} E. Lohrmann, M. W. Teucher, |! and 
Marcel Schein, ** Department of Physics, Uni- 
versity of Chicago, Chicago, Illinois (Received 
November 28, 1960). 


Eighty-four interactions of protons and neu- 
trons were located in a 22-liter stack of nuclear 
emulsion by tracing back showers of minimum- 
ionizing particles to their origins. The distri- 
bution of the number of shower particles, and 
the number of heavily ionizing prongs are pre- 
sented for 57 events with dip angles <17°. The 
average energy of these events is 3.510" ev. 
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The average number of shower particles emitted 
in nucleon-nucleon collisions at this energy is 
15+5, as estimated from 8 interactions without 
heavy prongs. The angular distributions of the 
shower particles are presented for the 57 events. 
They can be transformed into a system in 
which the angular distribution is roughly sym- 
metric. This is true even for the collisions with 
heavy target nuclei (Nj, >5). The degree of ani- 
sotropy of the angular distributions is in dis- 
agreement with a hydrodynamical model of 
nucleon-nucleus collisions. A lower limit for 
the collision mean free path of the primary par- 
ticles of 20 cm in emulsion was obtained. By 
scanning the forward cone of the primary inter- 
actions, 76 secondary interactions of charged 
and neutral shower particles were found. The 
distribution of the prong numbers, of the energy, 
and the characteristics of their angular distri- 
bution are presented. The best estimate of the 
ratio of secondary collisions produced by neutral 
particles, and the number produced by charged 
particles, is Ny/Ncep=0.4040.11. Adding this 
result to other published data, it is concluded 
that 30+ 6% of the particles produced in colli- 
sions having a primary energy of several Tev 
are not 7 mesons. A collision mean free path of 
41+8 cm was found for the forward-cone shower 
particles. 


*On leave from the Marquette University, Milwaukee, 
Wisconsin. 

TNow with Atomic Energy Commission, Bombay, 
India. 

ton leave of absence from the University of Buffalo, 
Buffalo, New York. 

Present address: Physikalisches Staatsinstitut, 

University of Hamburg, Hamburg, Germany. 

*Deceased. 


ANGULAR DISTRIBUTIONS OF SECONDARY 
PARTICLES PRODUCED IN HIGH-ENERGY 
NUCLEAR COLLISIONS AND THE TWO-CENTER 
MODEL OF MULTIPLE MESON PRODUCTION. 
J. Gierula,* D. M. Haskin, and E. Lohrmann, 
Department of Physics, University of Chicago, 
Chicago, Illinois (Received November 28, 1960). 


The angular distributions of shower particles 
from 54 nuclear interactions of protons and 
neutrons with energies > 10'* ev in a stack of 
nuclear emulsions are analyzed. The method 
consists essentially in normalizing the angular 
distributions of all events in the x =logtané@ 
scale to the same dispersion o. One finds a 
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very significant deviation from the normal dis- 
tribution predicted by hydrodynamical models. 
The deviation goes in the direction indicated by 
the two-center model (two maxima in the plot 
over the x coordinate). The correlation between 
the separation of the two emitting centers and o 
is also in qualitative agreement with the model. 
The angular distribution in the rest system of 
the emitting centers is found to be roughly iso- 
tropic. The two-center model also offers an 
explanation for certain characteristic features 
observed for the angular distribution of events 
with a small number of shower particles (n, <5). 
On the basis of this model a coefficient of in- 
elasticity of = 0.2 is obtained for these events. 

Interactions characterized by small evapora- 
tion (Nj, <5) and small numbers of shower parti- 
cles (n, < 20) show the characteristic two-maxi- 
mum shape. The same shape is found for, pre- 
sumably central, collisions with heavy nuclei 
in the emulsion (N;,>8, n,>40). However, the 
group of collisions with N;, <5, n,>20 does not 
show the two maxima. The last two observations 
cannot be explained by the present simple form 
of the two-center model. 

The results of this paper are in good agree- 
ment with a similar analysis carried out by 
Gierula, Miesowicz, and Zielinski. 


*On leave of absence from Cosmic Ray Department, 
Institute of Nuclear Research, Krakow, Poland. 


PHENOMENOLOGICAL STUDY OF PION PHOTO- 
PRODUCTION WITH POLARIZED PHOTONS. 

G. T. Hoff, Enrico Fermi Institute for Nuclear 
Studies, University of Chicago, Chicago, Illinois 
(Received December 6, 1960). 


Expressions for the angular distributions and 
polarizations of pion photoproduction with polar- 
ized photon beams are derived from a phenomeno- 
logical production matrix. The experiments 
necessary for the complete determination of the 
multipole amplitudes are discussed in general, 
and in particular for the case in which only con- 
tributions up to p waves in the final state are 
important. Complete determination requires 
circularly polarized beams. But if only s and p 
waves contribute, experiments with linearly 
polarized beams completely determine the pro- 
duction matrix. The knowledge of these ampli- 
tudes would allow the determination of the un- 
known p-wave scattering phase-shifts (a@,,, Q5, 
and a,,) up to energies of about 300 Mev. 
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Invariance properties of the angular distribu- 
tions and polarizations are found and tables are 
given. 


NUCLEAR CAPTURE OF MUONS WITH ELEC- 
TRON EMISSION. M. Conversi, L. di Lella, 

A. Egidi, C. Rubbia,* and M. Toller, Istituto 

de Fisica dell’ Universita, Sezione di Roma 
dell’ I. N. F. N., Roma, Italy (Received Decem- 
ber 15, 1960). 


Two experiments carried out to search for 
the process of muon capture with electron emis- 
sion are reported. The second of the two ex- 
periments is nearly 200 times more sensitive 
than earlier attempts to find this capture mode, 


but no indication is obtained in favor of the latter. 


In both experiments negative muons are made to 
stop in copper, where coherent capture is pre- 
dominant, so that the “capture electrons” should 
be emitted with an energy spectrum sharply 
peaked around 100 Mev. For the branching ratio 
of the process searched for, relative to ordinary 
muon capture, upper limits of about 5x107° and 
5x10~° are established, respectively, through 
the first and second experiment. 


*Now at CERN, Geneva, Switzerland. 


ELASTIC SCATTERING OF NEGATIVE PIONS 
BY PROTONS AT 230, 290, 370, AND 427 Mev. 
Lester K. Goodwin, Robert W. Kenney, and 
Victor Perez-Mendez, Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received December 5, 1960). 


The elastic differential cross section for the 
scattering of negative pions by hydrogen was 
measured at laboratory-system pion kinetic 
energies of 230, 290, 370, and 427 Mev. The 
elastically scattered pions were detected by a 
counter telescope which discriminated against 
recoil protons and inelastic pions on the basis 
of range. Differential cross sections were ob- 
tained at nine angles for each energy and were 
fitted by a least-squares program to a series of 
Legendre polynomials. At the three higher ener- 
gies, D waves are required to give satisfactory 
fits to the data. The real parts of the forward- 
scattering amplitudes calculated from this ex- 
periment are in agreement with the predictions 
of dispersion theory. The results of this ex- 


periment, in conjunction with data from other 
pion-nucleon scattering experiments, support 
the hypothesis of charge independence at these 
higher energies. 


PROPERTIES OF COSMIC-RAY AIR SHOWERS 
AT SEA LEVEL. G. W. Clark, J. Earl, W. L. 
Kraushaar, J. D. Linsley, B. B. Rossi, F. 
Scherb, and D. W. Scott, Laboratory for Nuclear 
Science, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received December 
13, 1960). 


An investigation at sea level of cosmic-ray air 
showers with sizes from 510° to over 10° par- 
ticles is described. The core locations, arrival 
directions, and particle density distributions of 
several thousand showers whose cores landed 
within an area of 10° m* were determined by the 
techniques of fast-timing and density sampling. 
The most important results are as follows: (1) 
The existence of primary particles with energies 
greater than 10'* ev is established by the obser- 
vation of one shower with more than 10° parti- 
cles; (2) the function 


f(r) = 0.45(N/R,2) r-°7(1+7)7*?, 


where r=R/R, and R,=79 m, describes the 
lateral distribution of particles at distances in 
the range 50 m<R<400 m and for showers with 
sizes in the range 5x10°< N<10®; (3) at dis- 
tances greater than 50 m from the core the den- 
sity fluctuations in individual showers have a 
Poisson distribution; (4) the size and zenith- 
angle distribution can be represented by the for- 
mula 


s(N,x) =s,(10°/ny" exp[- (x -x)/A], 


where x =x, cosé, x,=1040 g cm™, s,=(6.6+ 1.0) 
x10-® cm~? sec™! sr-*, F=1.940.1, A=(113+ 9) 
gem’, x,<x<1.3x,, and 7x10°<N<17x10*; (5) 
no evidence is found of anisotropy in the arrival 
directions or of a break in the energy spectrum 
of the primaries up to the largest energies ob- 
served; (6) assuming a specific model for show- 
er development and taking into account fluctua- 
tions in the depth of the first interaction, the in- 
tegral energy spectrum of the primaries is 

J(E) =J,(10"5/E)”, where J, =(8.1+ 3.1) x10-™ 
cm~ sec”! sr~', y=2.17, and 3x10"* ev<E< 10" 
ev. 
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RANGE OF PROTON-ANTIPROTON ANNIHILA- 
TION NEAR 1.0 Bev. Osamu Hara,* School of 
Physics, University of Minnesota, Minneapolis, 
Minnesota (Received December 9, 1960). 


The range of the proton-antiproton annihilation 
was calculated for antiprotons with energy near 
1 Bev. The point is to get the range of pure 
annihilation interaction, separating the effect 
of pion production. It was found that the root 
mean square of this range is given by (1.19+ 0.07) 
x107** cm almost independently of the energy. 


*On leave of absence from Nihon University, Tokyo, 
Japan. 


A BINDING IN HYPERNUCLEI BY NONLOCAL 
INTERACTION. G. Rajasekaran and S. N. Bis- 
was, Tata Institute of Fundamental Research, 
Bombay, India (Received October 21, 1960; 
revised manuscript received December 27, 1960). 


The characteristics of the A-N interaction at 
low energy have been obtained assuming that the 
A-N potential is nonlocal but separable and 
similar to that suggested by Yamaguchi in the 
case of an N-N potential. The unknown param- 
eters entering in the proposed potential are 
determined on the basis of the global symmetry 
hypothesis of the strongly interacting particles. 
Our model predicts, in agreement with Dalitz 
and Downs’ phenomenological findings of the 
nature of the A-N potential based on hyperfrag- 
ment data, that there is no bound A-nucleon 
system and that the singlet A-N potential is 
stronger than the triplet potential, both being 
attractive. Binding energies of the A particle 
in light hypernuclei based on the present model 
are, however, much too high compared with the 


experimental data. It is further pointed out that 
although the global symmetry hypothesis (g,, 
=£D_ =ZNz) Supplemented by the Yamaguchi- 
type nonlocal A-N potential is incompatible with 
the presently existing data, the restricted sym- 
metry (g,, =>, *8N_) model is certainly ad- 
missible. 


V - @ COLLISIONS IN THE LEE MODEL. R. D. 
Amado, Department of Physics, University of 
Pennsylvania, Philadelphia, Pennsylvania (Re- 
ceived November 16, 1960). 


The methods of dispersion theory are used to 
obtain an exact expression for the V - 6 elastic 
scattering amplitude and the amplitude for the 
production process, V+0~-N+26, in the Lee 
model. 


MAGNETIC MOMENTS OF THE A AND © HY- 
PERONS. Katsumi Tanaka, .Argonne National 
Laboratory, Argonne, Illinois (Received April 
25, 1960; revised manuscript received September 
20, 1960). 


A relation among the magnetic moments of ae 
=°, =, and A is obtained as a consequence of 
the proposed symmetries of strong interactions, 
a minimal electromagnetic coupling for the 
electromagnetic interactions being assumed. 
The magnetic moments of the A and = hyperons 
are calculated with the aid of mass spectral re- 
presentations in which only the contributions of 
the bound states are taken into account. The 
present calculations of these magnetic moments 
are compared with various other calculations. 
Remarks are made on the possible experimental 
values. 





